Interaction  between  Lung  Mechanics  and  Gas  Exchange  by  Low  Volume  Kig 
Frequency  Pulmonary  Ventilation  in  Patients  with  Respiratory  Failure 


Annual  Summary  Report 


Jeffrey  M.  Drazen,  M.D.,  Phillip  Drinker,  Ph.D.,  Thomas  Rossing,  M.D 
Julian  Solway,  M.D.,  Roger  D.  Kamm,  Ph.D.,  Isnmail  Dreshai,  M.D. 

and  Rayhana  Akhaven 


November  30,  1985 


Supported  by 


U.S.  ARMY  MEDICAL  RESEARCH  AND  DEVELOPMENT  COMMAND 
Fort  Detrick,  Frederick,  Maryland  21701-5012 


Contract  No.  DAMD17-82-C-2210 


Brigham  and  Women's  Hospital 
75  Francis  Street 
Boston,  Massachusetts  02115 


DOD  Distribution  Statement 


Approved  for  public  release;  d istribuiion  unlimited 


ft*4. 


fedllftlTY  TLAjSlFI&ATION  OP  THIS  PAC 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassif ied 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 


0<62 


REPORT  DOCUMENTATION  PAGE 

lib  RESTRICTIVE  MARKINGS 


2b.  OECLASSIFICATION  /  DOWNGRADING  SCHEDULE 
T  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL 
Brigham  &  Women's  Hospital  (if  applicable) 

6c.  ADDRESS  (Cty,  State,  and  IIP  Code) 

10  Vining  Street 
Boston,  MA  02115 

8a.  NAME  OF  FUNDING  /  SPONSORING  8b.  OFFICE  SYMBOL 

ORGANIZATION  U.S.  Army  Medical  (>f  applicable) 
Research  &  Development  Command 
8c  ADDRESS  (City,  State,  and  ZIP  Code) 

Fort  Detrick 

Frederick,  Maryland  21701-5012 


:  3  DISTRIBUTION/ AVAILABILITY  OF  REPORT 


fS.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 


|  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

DAMD17-82-C-2210 


10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  PROJECT 

ELEMENT'  NO.  NO. 

62734A'  3M162734A87 


PROJECT- 

TASK, 

NO 

NO. 

3M162734A87 

5  CC 

WORK  UNIT 
ACCESSION  NO. 


11  TITLE  (Include  Security  Classification) 

Interaction  Between  Lung  Mechanics  and  Gas  Exchange  by  Low  Volume  High  Frequency 
Pulmonary  Ventllat-irm  in  Pat.ipnts  with  BocptT-ai-rM-y  FaiTur^  _  . -  -  - - 

12  PERSONAL  AUTHOR(S) 

Jeffrey  Drazen.  M„D . 

13a.  TYPE  OF  REPORT  13b.  TIME  COVERED  14.  DATE  OF  REPORT  (Year,  Month,  Day)  IS  PAGE  COUNT 

Annual _  from  10/1/84  Toq/WSSl  85  November  70 

16  SUPPLEMENTARY  NOTATION 

17  _ COSATi  COOES _  18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  Nock  number) 

FIELD _ GROUP _ SUB-GROUP  High  Frequency  Ventilation 

__15 _ 06 _  _ I  Mechanical  Ventilation 

06  r.  n  _.  l  .  _  i - 

19  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  Nock  number) 

Flow  limitation  during  a  forced  expiration  was  simulated  by  a  mathematical 
model.  This  model  draws  upon  a  pressure-area  law  obtained  in  other  work,  and 
ou  known  methods  of  analysis  for  flow  in  collapsible  tubes.  This  approach 
represents  an  improvement  over  previous  models  in  that  (i)  the  effects  of 
changing  lung  volume  and  of  parenchymal/ bronchial  interdependence  are 
simulated,  (ii)  a  more  realistic  representation  of  the  collapsed  airways  is 
employed,  (iii)  a  solution  is  obtained  mouthward  of  the  flow  limiting  site  by 
allowing  for  a  smooth  transition  from  subcritical  to  supercritical  flow  speeds, 
then  matching  mouth  pressure  by  imposing  an  elastic  jump  (an  abrupt  transition 
from  supercritical  to  subcritical  flow  speeds)  at  the  appropriate  location,  and 
(iv)  the  effects  of  levels  of  effort  (or  vacuum  pressure)  in  excess  of  those 
required  to  produce  incipient  flow  limitation  are  examined,  including  the 
effects  of  potential  physiological  limitation. 

20.  DISTRIBUTION  /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  DTIC  USERS 
22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 

MrsoVir^T^^^Hcr^^ ^ 

DO  FORM  1473,  84  MAR  83  APR  edition  may  be  used  unt.l  exhausted.  SECURITY  CLASSIFICATION  OF  THIS  PA; 

All  other  editions  are  obsolete. 


□  DTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

22b  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL  ““ 
til  exhausted.  SECURITY  CLASSIFICATION  OF  THIS  PAGE 


SECURITY  CLASSIFICATION  OF  THIS  RAGE 


In  patients,  dynamic  hyperinflation  of  the  lungs  occurs  during  high  frequency 
oscillatory  ventilation  (HFOV),  and  has  been  attributed  to  asymmetry  of 
inspiratory  and  expiratory  impedances.  To  identify  the  nature  of  this 
asymmetry,  we  compared  changes  in  lung  volume  (Vl)  observed  during  HFOV  in 
ventilator-dependent  patients  with  predictions  of  Vl  changes  from  electrical 
analogs  of  3  potential  modes  of  impedance  asymmetry.  In  the  patients  when  a 
fixed  oscillatory  tidal  volume  was  applied  at  a  low  mean  airway  opening 
pressure  (Pao)  which  resulted  in  little  increase  in  FRC,  progressively  greater 
dynamic  hyperinflation  was  observed  as  HFOV  frequency  ( f )  was  increased.  When 
mean  Pao  was  raised  such  that  resting  Vl  increased,  Vl  remained  at  this 
level  during  HFOV  as  f  was  increased,  until  a  critical  f  was  reached;  above 
this  value,  Vl  increased  further  with  f  in  a  fashion  nearly  parallel  to  that 
observed  when  low  mean  Pao  was  used.  Three  modes  of  asymmetric  inspiratory  and 
expiratory  impedance  were  modelled  as  electrical  circuits:  (a)  fixed 
asymmetric  resistance  [Rexp  >  Rinsp];  (b)  variable  asymmetric  resistance 
[Rexp(VL)  >  Rinsp,  with  Rexp(VL)  decreasing  as  Vl  increased;  and  (c) 
equal  Rinsp  and  Rexp,l  but  with  superimposed  expiratory  flow  limitation,  the 
latter  simulated  using  a  bipolar  transistor  as  a  descriptive  model  of  this 
phenomenon.  The  fixed  and  the  variable  asymmetric  resistance  models  displayed 
a  progressive  increase  of  mean  Vl  with  f  at  either  low  or  high  mean  Pao. 

Only  the  expiratory  flow  limitation  model  displayed  a  dependence  of  dynamic 
hyperinflation  on  mean  Pao  and  f  similar  to  that  observed  in  our  patients.  We 
conclude  that  expiratory  flow  limitation  can  account  for  dynamic  pulmonary 
hyperinflation  during  HFOV. 
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I.  SUMMARY  j 

Flow  limitation ^tiring  a, forced  expiration  was  simulated i by  a  mathematical 
model.  This  model  draws  upon  a  pressure-area  law  obtained  in  other  work,  and 
on  known  methods  of  analysis  for  flow  in  collapsible  tubes.  This  approach 


represents  an  /improvement  over  previous  models  in  that 


the  effects  of 


changing  lung/ volume  and  of  parenchymal/ bronchial  interdependence  are 
simulated,  (~±±)/a  more  realistic  representation  of  the  collapsed  airways  is 
employed,  a  solution  is  obtained  mouthward  of  the  flow  limiting  site  by 

allowing  for  a  smooth  transition  from  subcritical  to  supercritical  flow  speeds, 
then  matching. mouth  pressure  by  imposing  an  elastic  jump  (an  abrupt  transition 
from  supercritical  to  subcritical  flow  speeds)  at  the  appropriate  location,  and 
the  effects  of  levels  of  effort  (or  vacuum  pressure)  in  excess  of  those 
required  to  produce  incipient  flow  limitation  are  examined,  including  the 
effects  of  potential  physiological  limitation.  ^  k 

In  patients,  dynamic  hyperinflation  of  the  lungs  occurs'7 during  high  frequency 
oscillatory  ventilation  (HFOV),  and  has  been  attributed  to  asymmetry  of 
inspiratory  and  expiratory  impedances.  To  identify  the  nature  of  this 
asymmetry,  we  compared  changes  in  lung  volume  (V£)  observed  during  HFOV  in 
ventilator-dependent  patients  with  predictions  of  V£  changes  from  electrical 
analogs  of  3  potential  modes  of  impedance  asymmetry.*  In  the  patients  when  a 
fixed  oscillatory  tidal  volume  was  applied  at  a  low  mean  airway  opening 
pressure  (Pao)  which  resulted  in  little  increase  in  FRC,  ^progressively  greater 
dynamic  hyperinflation  was  observed  as  HFOV  frequency  ( f )  Was  increased.  When 
mean  Pao  was  raised  such  that  resting  Vl  increased,  VL  remaided  at  this 
level  during  HFOV  as  f  was  increased,  until  a  critical  f  was  reached;  above 
this  value,  Vl  increased  further  with  f  in  a  fashion  nearly  parallel  to  that 
observed  when  low  mean  Pao  was  used.  Three  modes  of  asymmetric  inspiratory  and 
expiratory  impedance  were  modelled  as  electrical  circuits:  (a)  fixed 
asymmetric  resistance  [Rexp  >  Rinsp];  (b)  variable  asymmetric  resistance 
[Rexp(VL)  >  Rinsp,  with  Rexp(VL)  decreasing  as  Vl  increased;  and  (c) 
equal  Rinsp  and  Rexp,l  but  with  superimposed  expiratory  flow  limitation,  the 
latter  simulated  using  a  bipolar  transistor  as  a  descriptive  model  of  this 
phenomenon.  The  fixed  and  the  variable  asymmetric  resistance  models  displayed 
a  progressive  increase  of  mean  Vl  with  f  at  either  low  or  high  mean  Pao. 

Only  the  expiratory  flow  limitation  model  displayed  a  dependence  of  dynamic 
hyperinflation  on  mean  Pao  and  f  similar  to  that  observed  in  our  patients.  We 
conclude  that  expiratory  flow  limitation  can  account  for  dynamic  pulmonary 
hyperinflation  during  HFOV. 
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III.  OVERALL  INTRODUCTION 


t 

During  year  03  of  this  contract  we  performed  both  hardware  model  j 

experiments  and  patient  studies.  Each  will  be  addressed  separately. 


1.  INTRODUCTION 

A  forced  expiration  maneuver,  as  characterized  by  a 
Maximal  Expiratory  Flow  Volume  (MEFV)  curve  has  long  been  used 
as  a  test  of  pulmonary  function.  (See  Hyatt  (1983)  for  an 
extensive  review.)  The  MEFV  curves  produced  by  a  given  subject 
are  highly  reproducible,  but  can  vary  markedly  among  different 
subjects  even  when  an  adjustment  is  made  to  correct  for 
differences  in  lung  volume  (Green  et.  al. ,  1974).  Consistent 
and  characteristic  differences  in  the  shape  of  the  MEFV  curve 
are  found  among  different  species  (Macklem  et.  al. ,  1968)  and 
depending  on  whether  the  forced  expiration  is  produced  by 
respiratory  effort  or  by  lowering  mouth  pressure  (reference). 

The  repeatibility  of  the  MEFV  curve  in  a  particular 
subject  is  often  attributed  to  flow  limitation;  that  the  flow 
rate  produced  at  constant  lung  volume  is  insensitive  to  further 
changes  in  respiratory  effort  when  the  former  reaches  some 
critical  level.  For  lung  volumes  from  about  30%  to  70%  of 
Total  Lung  Capacity  (TLC)  this  is  generally  thought  to  be  due 
to  wave  speed  limitation  which  occurs  when  the  flow  velocity 
equals  the  wave  speed  at  some  location  in  the  airways.  The 
point  at  which  these  speeds  are  equal  (termed  the  Flow  Limiting 
Site  (FLS) )  is  apparently  found  in  the  trachea  at  high  lung 
volumes,  then  progresses  upstream  (toward  the  periphery)  as 
lung  volume  decreases  (Smaldone  et  al.,  1976;  Hyatt  etal., 

1980;  Pedersen  et  al.,  1982),  but  probably  resides  in  the 
central  airways  even  at  low  lung  volumes  (Smaldone  and  Smith, 
1985).  While  these  arguments  seem  essentially  correct  for  most 
lung  volumes,  at  lung  volumes  greater  than  about  70%  TLC, 
maximum  flow  is  determined  instead  by  physiological  limitation 
of  the  expiratory  muscles  (Hyatt  et  al.,  1966;  Mead  et  al., 
1967)  in  which  case  expiratory  flow  is  observed  to  be  effort 
dependent  (Mead  et  al,1967). 

Recently,  experimental  methods  have  been  developed  to 
determine  the  location  of  the  FLS  in  excised  and  intact  canine 
and  human  lungs.  In  these  studies,  the  FLS  is  taken  to  be  the 
furthest  mouthward  point  at  which  lateral  pressure,  at  a 
certain  lung  volume,  remains  constant  while  mouth  pressure  is 
lowered  in  successive  expirations.  Most  evidence  suggests  that 
the  FLS  first  forms  in  the  trachea,  then  progresses  upstream 
(peripherally),  with  decreasing  lung  volume. 

The  nature  of  flow  limitation  has  also  been  examined 
through  experiments  using  different  gas  mixtures.  The  MEFV 
curves  are  found  to  depend  primarily  upon  gas  density  at  high 
lung  volumes,  but  become  Increasingly  sensitive  to  gas 
viscosity  as  lung  volume  decreases.  The  relative  effects  of 
different  gases  on  the  location  of  the  FLS  has  not  been  clearly 
established  (Wood  et  al.,  1976;  Mink  et  al. ,  1979;  Staats  et 
al.,  1980). 
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These  observations,  and  the  desire  to  better  understand 
the  mechanistic  basis  for  flow  limitation,  provided  the  impetus 
to  develop  models  to  simulate  a  forced  expiration.  The  most 
recent  and  most  rigorous  of  these  models  is  that  of  Lambert  and 
Wilson  .  In  their  model,  the  lung  is  treated  as  a  symmetric 
branching  network  with  distributed  properties.  Flow 
calculations  proceed  from  the  alveolar  zone  toward  the  mouth 
using  differential  forms  of  the  mass  and  momentum  conservation 
equations  in  conjunction  with  a  constitutive  relation  between 
tranpulmonary  pressure  and  cross-sectional  area  (the  "tube 
law").  Flow  limitation  is  attained  by  incrementing  the  flow  at 
the  upstream  (alveolar)  end  of  the  system  until  reaching  a 
state  is  reached  for  which  the  flow  speed  and  wave  speed  become 
equal  at  some  point  in  the  network  (the  FLS).  While  their 
model  reproduces  many  features  of  a  forced  expiration,  it 
cannot  be  used  to  examine  flows  downstream  of  the  FLS. 

Further,  the  airway  constitutive  relation  used  in  their  study 
is  based  on  limited  data  obtained  by  uniform  lung  inflation  at 
positive  distending  pressures  and  fails  to  take  into  account 
the  effects  of  parenchyma/bronchial  interdependence  or  changes 
in  lung  volume. 

A  more  general  approach  to  steady  flow  in  collapsible 
tubes  was  developed  by  Shapiro  (1977).  Starting  with  the 
equations  of  mass  and  momentum  conservation  and  the  tube  law 
for  the  vessel,  the  analysis  shows  how  the  area  ratio  (a  =  A/A 
where  A  is  some  reference  cross-sectional  area)  and  the  speed0 
index  (5,  the  ratio  between  flow  velocity  and  the  local  speed 
of  long  waves)  vary  with  streamwlse  distance.  This  formulation 
more  clearly  identifies  the  different  characters  of  sub- 
critical  (S< 1)  and  super-critical  (S>1)  flow  and  allows  one  to 
examine  possible  transitions  between  these  two  states. 

In  a  previous  work  (Elad  et  al. ,  1985)  we  employed  this 
same  approach  to  investigate  various  types  of  flow  in  a  lung¬ 
like  model.  Solutions  were  obtained  which  contained  a  smooth 
sub-  to  super-critical  transition,  followed  by  a  "shock-like" 
deceleration  (an  "elastic  jump")  back  to  sub-critical  speed. 

The  location  and  strength  of  the  elastic  jump  provided 
solutions  which  matched  "mouth"  pressure  downstream.  These 
studies,  while  qualitatively  similar  to  flows  during  a  forced 
expiration,  were  conducted  for  the  primary  purpose  of  exploring 
different  types  of  flow  behavior.  Accordingly,  no  attempt  was 
made  to  mimic  all  important  aspects  of  the  respiratory  system. 

In  the  present  investigation  we  employ  the  same  principles 
to  a  model  which  more  closely  simulates  the  lung.  This  model 
allows  us  to  calculate  flow  conditions  beyond  the  FLS,  an 
essential  step  toward  the  development  of  models  incorporating 
airway  asymmetry.  We  also  consider  the  conditions  that  exist 
when  respiratory  effort  exceeds  the  minimum  value  necessary  to 
produce  flow  limitation  and  demonstrate  how  the  occurance  of 
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"elastic  jumps"  allows  the  flow  solution  to  match  the 
downstream  pressure  existing  at  the  mouth. 

These  calculations  incorporate  a  tube  law  based  on 
experimental  results  in  which  bronchial  pressure  is  varied 
while  alveolar  pressure  is  held  constant.  This  tube  law 
provides  for  a  more  realistio  representation  of  the  effects  of 
parenchymal  tensions  and  of  changing  lung  volume  as  the  forced 
expiration  proceeds  to  completion. 

The  differences  and  similarities  between  forced  expiration 
and  forced  evacuation  are  examined  as  are  the  effects  of  small 
perturbations  from  a  monotonic  variation  in  airway  stiffness  or 
area,  as  might  exist  in  various  respiratory  disorders. 


In  Elad  et  al  (1985),  a  generalized,  self -similar  tube  law 
was  obtained,  based  on  the  premise  that  the  fundamental  character 
of  collapse  at  different  lung  volumes  is  essentially  the  same. 
Accordingly,  we  found  that  an  appropriately  chosen  dimensionless 
form  of  the  pressure-area  relationship  remains  valid  over  a  broad 
range  of  lung  volumes.  This  relationship  was  expressed  in  the 
form: 


(pbr-paxv>-p0 


a0-5  -  a'0-2 


where 


Here,  P,  is  the  internal  bronchial  pressure,  P  ,  the  alveolar 
pressure?  A  the  cross-sectional  area  of  a  singleairway ,  AQ  and 
PQ  are  the  cross-sectional  area  and  the  pressure  difference  (P, 
-P  i  )  at  the  inflection  point  of  the  pressure-area  curve,  and  K 
isatlie  effective  wall  stiffness.  In  Elad  et  al  (1985),  Aq,  Pq, 
and  K  were  shown  to  depend  upon  a  dimensionless  lung  volume,  >(, 
defined  by 


A  -  VVLO 


where  V.  is  lung  volume  and  V_  n  a  reference  volume  corresponding 
to  35%  Cf  TLC  (Klingele  et  alT ;  1971;  Hughes  et  al.,  1972). 
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We  now  need  to  generalize  from  the  tube  law  of  Elad  et  al 
(1985)  (which  was  based  on  data  from  a  single  airway  generation) 
to  one  which  is  applicable  to  the  entire  airway  system.  Again  we 
start  with  the  assumption  that  the  mechanics  of  airway  collapse 
is  basically  the  same  regardless  of  airway  size  or  location  in 
the  lung.  Therefore,  we  anticipate  the  existence  of  a  universal 
tube  law,  of  the  same  basic  form  of  equation  (1),  but  in  which 
the  parameters  which  characterize  the  stiffness  and  size  of  the 
airway  vary  as  a  function  of  airway  generation.  Thus  we  take  K, 

P  ,  and  A  to  be  functions  of  both  lung  volume  and  distance  along 
tne  branching  network.  Assuming  the  dependence  on  lung  volume  to 
be  of  the  same  form  for  all  airways,  we  carry  over  the  results  of 
Elad  et  al  (1985): 


p0  =  Poo(°  C0-0136  7,4  -  (4) 

Aq  -  AQ0(O  Cl.O  +  0.01345  (A2'7  -  1)]  (5) 

K  =  KQ(0  [  1.0  +  0.078  (A2,7  -  1)]  (6) 


where  pq0(£)  and  Aoo^^  ar8,  z'Qspsotiveiy .  the  distributions  of 
pressure  and  cross-sectional  area  at  the  inflection  point  of  the 
pressure-area  tube  law  for  lung  volume  VTn,  K  (£)  is  similarly 
the  distribution  of  effective  wall  stiffness  at  V_Q,  and  £  -  x/L 
is  a  nondimensional  length,  normalized  by  the  total  length  L  of 
the  conducting  system  up  to  the  17th  generation. 


Changes  Along  the  Airway  Network 

According  to  the  discussion  of  Elad  et  al  (1985),  AQ 
represents  the  distribution  of  area  at  the  Inflection  point  of 
the  tube  law  for  that  lung  volume  at  which  the  parenchyma  exerts 
no  external  tethering  force,  a  value  assumed  earlier  to  be  equal 
to  Vj_.  In  order  to  use  the  morphological  data  of  Weibel  (1963), 
whicn  was  obtained  from  a  lung  inflated  to  75%  TLC ,  we  need  to 
adjust  down  to  V-0  (35%  TLC),  taking  into  account  the  associated 
changes  in  airwajrlength.  For  this  purpose,  we  assume  that  the 
length  of  each  airway  is  constant  for  changes  at  constant  lung 
volume,  but  varies  as  the  cube  root  of  the  volume  ratio  (Hughes 
et  al. ,  1972).  Hence 


W 


1/3 


(7) 


where  LQ0  is  the  length  of  the  vessel  at  VLQ.  Taking  into 
accountboth  the  change  in  lung  volume  and^the  length  correction 
gives 


A00CO 


0.916  Aw(£) 


(8) 


where  A^(£)  is  the  area  distribution  of  Weibel  (1963)  given  later 
in  equation  (9). 
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It  is  somewhat  more  difficult  to  obtain  accurate  estimates 
for  the  two  parameters  affecting  airway  stiffness,  K  and  Pqq. 

Many  observations  confirm  that  the  peripheral  airways  are  most 
compliant  and  that  the  airways  become  progressively  stiffer 
moving  toward  the  trachea.  There  are,  however,  no  data  available 
from  which  we  can  determine  the  precise  nature  of  this 
distribution.  Hence,  we  assume  that  the  effective  wall  stiffness 
at  VLQ,  varies  exponentially  along  the  bronchial  tree,  according 
to 


Ko(£)  -  K00  e*P(Ck£) 


cm 


where  is  the  effective  wall  stiffness  of  the  most  peripheral 
generation  (generation  17  in  the  present  model)  at  VLfl  ,  and  C. 
is  an  empirical  constant  that  controls  the  relative  change  along 
the  bronchial  generations. 

We  obtain  one  value  for  K  from  the  data  analysis  of  Elad  et 
al  (1985),  and  set  K( £=0 . 33 , A-i)  -  37.1  cmH„0  (data  from  the 
third  generation  of  Dog  20).  Since  no  data^of  comparable  form 
exist  for  airways  from  other  generations,  we  consider  instead  the 
experimental  results  of  Martin  et  al.,  (1958),  obtained  from 
three  excised  bronchi:  trachea,  5mm  and  2mm  outside  diameter. 
Noting  that  the  2mm  airway  is  roughly  6  times  as  compliant  as  the 
trachea  we  obtain  C. -2.40  which,  in  turn,  dictates  that  Knn  - 
16.7cmH  0.  *  uu 

Finally,  we  need  to  specify  the  variation  in  PQ.  Here,  for 
lack  of  empirical  data  on  which  to  rely,  we  assume  that  P-/K  is 
independent  of  £  (l.e. ,  PQ  varies  exponentially)  for  reasons  of 
geometric  and  structural  similarity. 

Airway  System  Geometry 

For  the  purpose  of  the  present  model,  the  bronchial  tree  is 
represented  as  a  continuously  branching  symmetric  network  with  a 
monotonic  increase  of  the  total  cross-sectional  area  (Fig.l) 
based  on  the  morphometric  data  of  Weibel  (1963),  obtained  from  a 
human  lung  Inflated  to  75%  of  TLC.  The  complex  geometry  with 
discrete  bifurcations  is  replaced  by  a  continous  one-dimensional 
model  in  which  the  number  of  the  bronchi  and  their  diameter  at 
each  generation  is  represented  by  smooth  and  continuous  functions 
of  distance  along  the  bronchial  tree. 

We  describe  the  total  cross-sectional  area  (A)  and  number  of 
branches  (N)  using  hyperbolic  functions  with  an  exponential 
correction  term,  fitted  to  Weibel 's  data.  The  result  is  shown  in 
Fig  2.,  and  is  given  by: 


Aw(0  -  2.3255(00.01)  0,9  -  14 . 0exp(-3 . 98440 
Ny(0  -  1.0387(^+0. Ol)-2,4  -  50. 0exp(-5. 97660 


(10) 
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Also,  from  Weibel,  we  find  that  L  =  26.5  cm  at  75%  TLC,  or  20.6 
cm  at  our  reference  state  according  to  equation  (2).  The  local 
diameter  of  a  single  tube  at  a  given  lung  volume  is  given  by 


Do 


7TNC0 


Consistant  with  other  investigations  Ce.g.  Lambert  and 
Wilson) ,  we  assume  that  flow  during  a  forced  expiration  can  be 
described  by  a  series  of  quasi-steady  states  at  progressively 
decreasing  fixed  lung  volumes.  This  is  justified  by  the  fact 
that  the  characteristic  time  for  a  forced  expriation  is  much 
larger  than  both  the  time  required  for  waves  to  propagate  along 
the  bronchial  tree  and  for  a  fluid  particle  to  traverse  the 
branching  network.  We  consider  a  one-dimensional  flow  with  wall 
friction  in  a  multi-channel  network.  External  pressure  is 
assumed  to  be  equal  to  alveolar  pressure,  as  described  in  Elad  et 
al  (1985),  We  later  examine  the  effect  of  a  gradient  in  external 
pressure  at  the  transition  between  the  intrapulmonary  and 
extrapulmonary  airways.  With  these  assumptions,  we  can  write  the 
equations  for  mass  conservation: 


uA  **  constant 


(12) 


and  conservation  of  momentum 


1  dPx 


(13) 


where  s  is  the  airway  perimeter,  t  the  shear  stress  at  the  wall, 
u  the  cross-sectional  average  gas  velocity,  A  the  cross-sectional 
area  of  the  tube,  p  the  gas  density  and  x  the  longitudinal 
coordinate . 


The  shear  stress  is  expressed  as  t  u  f  J2  where  f_  is  the 
friction  coefficient.  The  friction  coefficient  is  determined 
from  Reynolds  (1982),  who  measured  the  total  frictional  loss  due 
to  water  flow  in  a  latex  model  of  the  first  10  generations  of  a 
human  bronchial  tree.  Based  on  his  experimental  results,  and  an 
assumed  distribution  of  loss  among  the  generations,  he  proposed 
the  following  formula  for  the  local  frictional  pressure  drop: 


V  «-  .*  V  -’ 


(a  +  b  Re)  4P, 


(14) 


where  is  the  Poiseuille  pressure  drop  87r/iQL/A  ,  Re  is  the 
Reynolds^number  (R  -  uD  /i--  where  v  is  the  kinematic  viscosity 
experimental  data  ®a  -  1:9,  b  -  0.0035).  This  result  can  also  be 
expressed  in  terms  of  the  friction  coefficient: 


16 

—  (a  +  b  Re) 
Re 


(15) 


which  exhibits  laminar  flow  behavior  when  Re  is  small  and  takes 
on  the  character  of  turbulent  flow  when  Re  is  large. 

We  have  chosen  to  use  this  relation  in  a  slightly  modified 
form,  using  b  =  0.0013  rather  than  0.0035  as  given  by  Reynolds. 
The  value  given  by  Reynolds  for  b  is  roughly  four  times  that 
found  in  fully  developed  turbulent  flow  through  a  moderately 
rough  straight  tube.  Based  on  the  recognition  that  (i)  Re  at  the 
point  of  greatest  frictional  loss  is  of  order  10  ,  and  (il)  that 
the  airways  where  this  loss  occurs  are  likely  to  be  partially 
collapsed  and  therefore  of  a  different  cross-sectional  shape  from 
the  airways  tested  by  Reynolds,  we  use  the  lower  value  given 
above.  Later,  in  the  Discussion,  we  examine  the  effect  of 
varying  both  a  and  b  on  the  simulation  results. 


Normalized  relationships  can  be  obtained  by  introducing  the 
speed  index,  S-u/c,  the  area  ratio  a  defined  in  equation  (2),  a 
normalized  alveolar  pressure  77  ,  -P  ,  /K,  a-nd  the  non-dimensional 
length  £”X/L .  The  resulting  expressions  are  similar  to  those 
obtained  in  our  previous  work  (Elad  et  al.  1985b) 
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where  II'  =  d/7/da  and.  M  -  all' ' /II'  .  This  presentation  highlights 
the  separate  effects  associated,  respectively,  with  variations  in 
rest  area,  variations  in  external  (alveolar)  pressure,  friction, 
and  variations  in  wall  stiffness  represented  by  the  first  through 
fourth  terms  on  the  RHS  of  each  expression. 

It  can  be  seen  that  mathematical  singularities  appear  when 
S-l  since  the  denominator  on  the  RHS  of  both  expressions  goes  to 
zero.  It  has  been  shown  (Shapiro,  1977,  Elad  et  al,  1985b)  that 
one  of  two  situations  can  exist  when  S-l.  If  the  numerator 
passes  through  zero  at  precisely  the  same  point  at  which  S-l  then 
a  smooth  and  continuous  transition  from  subcritical  to 
supercritical  speed  is  possible.  This  can  occur  in  a  variety  of 
ways.  Consider,  for  example,  the  case  in  which  only  the  first 
term  on  the  RHS  of  equations  (16)  and  (17)  is  important,  that 
associated  with  variations  in  A  .  A  smooth  transition  in  this 
instance  can  result  only  when  A°  passes  through  a  local  minimum 
(dA  /df  -  0)  when  S-l.  Under°more  general  circumstances  such 
as  are  found  in  the  lung,  all  the  terms  on  the  RHS  contribute, 
but  their  sum  must  still  pass  through  zero  when  S-l  in  order  to 
allow  a  smooth  transition.  While  it  would  seem  to  be  difficult 
to  achieve  such  a  singular  condition,  it  proves  to  be  relatively 
easy  to  produce  smooth  sub~to  supercritical  transitions  in 
laboratory  experiments. 

If,  on  the  other  hand,  the  numerator  on  the  RHS  does  not 
pass  through  zero  when  S-l,  a  different  situation  exists,  termed 
a  "choke  point."  At  a  choke  point,  since  the  numerator  is  non¬ 
zero,  the  derivatives  dS  /d£  and  da/d£  are  infinite  in  magnitude 
when  S  reaches  1.  For  related  reasons,  it  is  not  possible  to 
continue  the  mathematical  solution  downstream,  beyond  a  choke 
point.  Hence,  in  a  real  flow,  a  choke  point  can  only  occur  at 
the  end  of  the  airway  system.  A  FLS  must  therefore  be  a  sub-  to 
super-critical  transition  point  rather  than  a  choke  point . 

This  raises  a  question  concerning  the  nature  of  flow 
downstream  of  the  transition  or  FLS.  Following  a  smooth 
transition  to  S>1,  the  flow  is  supercritical  for  some  distance. 

In  laboratory  experiments  ,  this  super-critical  zone  can  be 
terminated  in  an  'elastic  jump'  if  the  downstream  pressure  is  too 
high;  this  is  a  shock-like  transition  back  to  subcritical  speed 
associated  with  a  relatively  abrupt  increase  in  area  and 
pressure.  Elastic  jumps  have  been  the  subject  of  several 
modeling  efforts  which  have  arrived  at  the  following 
description.  Application  of  conservation  of  mass  and  momentum 
for  the  region  containing  the  elastic  jump,  produces  the  result 
that  the  quantity 
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has  the  same  value  at  both  ends  of  the  elastic  jump 


The  physics  of  a  forced  expiration  is  controlled  by  Eqs. 
(18)  and  (19)  which  include  the  fluid  flow  dynamics  (Eqs.  (16) 
and  (17))  and  the  tube  wall  properties  (Eqs. (6),  (7),  (8),  (9), 
(10),  (11))  and  can  be  integrated  numerically  for  given  boundary 
conditions.  In  the  case  of  a  forced  expiration  mouth  pressure 
(at  £-1)  equals  atmospheric  pressure.  At  the  periphery  (£-0), 
the  bronchi  merge  into  the  alveolar  zone.  Consequently  the 
boundary  conditions  are : 


P  =P 
*br  alv 


(19) 


where  P  .  is  atmospheric  pressure  and  A?  is  the  pressure  drop 
across  glottis  and  mouth.  ® 

To  determine  ^P  ,  we  use  the  results  of  Jaeger  et  al. 
(1969)  who  demonstrated  that  upper  airway  resistance,  from  mid¬ 
trachea  to  the  mouth,  can  be  described  over  the  Reynolds  number 
range  Re-10  to  10**  by  a  nondimensional  loss  coefficient: 
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We  fitted  Jaeger's  result  for  C.  to  an  analytic  function  of 
Reynolds  number  (shown  in  Fig.  3)  given  by 

InC .  - 3.394  +  0.776  ln„Re  -  0.0495  (In  Re)2  + 

0.000813  (In  Re) 


(21) 


where  we  have  assumed,  as  did  Jaeger  et  al  _  that  cross-sectional 
area  at  the  glottis  is  approximately  1.0  cm  . 

In  addition,  it  has  been  shown  that  the  resistance  of  the 
upper  airways,  including  glottis  and  mouth,  decreases  with 
increasing  lung  volume  (Hyatt  et  al. ,  1961;  Wilson  1969; 

Stanescu  et  al.,  1972).  For  the  present  wort  we  assume  that  the 
lung  volume  dependence  is  of  the  form  of  4P./Q-1.5/A  which  is 
roughly  consistent  with  experimental  observation.  The  final 
expression  for  the  pressure  drop  across  the  glottis  and  the  mouth 
becomes 
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As  noted  above  the  FLS  is  a  transition  point  corresponding 
to  S-l.  Setting  S  -  1  in  Eqs.  (16)  and  (17)  gives  rise  to  a 
mathematical  singularity  in  that  the  denominator  of  both 
expressions  goes  to  zero.  For  a  mathematical  solution  to  exist 
on  both  sides  of  this  singular  point,  it  is  necessary  that  the 
denominator  of  both  expressions  also  equal  zero  when  S-l; 
consequently  da/df-dS'VdS-O/O .  Only  under  these  conditions  can  a 
continuous  solution  exist  which  causes  the  flow  to  accelerate 
from  S<1  to  S>1.  One  can  calculate  all  possible  singular 
(transition)  points  by  setting  the  numerators  of  the  right-hand 
side  of  Eqs.  (16)  and  (17)  (which  are  identical  for  S-l)  to  zero 
while  S-l.  This  yields  a  locus  of  possible  transition  points, 
for  each  lung  volume,  as  shown  in  Fig.  4.  A  transition  curve  for 
each  lung  volume  and  £*  may  be  found  by  integrating  Eqs.  (16)  and 
(17)  from  the  critical  point  (£  )  to  the  upstream  end  (£-0).  The 
relevent  curve  for  a  particular  l^ng  volume,  found  by  trial  and 
error,  is  that  for  the  value  of  £  that  satifies  the  boundary 
condition  P.  -p  .  at  £-0.  These  values  of  £*  representing  the 
flow  conditionsaare  indicated  by  the  open  circles  (corresponding 
to  each  of  several  different  lung  volumes)  on  the  curves  of  Fig. 
4.  Having  identified  the  FLS,  one  can  easily  calculate  the 
critical  flow  rate  (maximum  flow  rate)  for  a  given  lung  volume. 
Then,  by  integration  of  Eqs.  (18)  and  (19),  both  upstream  and 
downstream  (curves  A  and  B,  respectively,  in  Fig.  xx),  we  can 
obtain  S(£)  and  <*(£)  for  each  lung  volume. 

Selection  of  the  transition  point  as  shown  above  insures 
that  the  upstream  boundary  condition  will  be  satisfied.  But  this 
procedure  generally  leads  to  some  pressure  other  than  mouth 
pressure  at  £=0.  To  satisfy  the  downstream  boundary  condition  of 
pressure  requires  a  trial  and  error  procedure.  First  we 
integrate  eqns.  (18)  and  (19)  to  get  S(£)  and  a(£)in  the 
upstream  direction,  starting  from  £-1  where  we  know  both  the  flow 
rate  and  the  local  (mouth)  pressure.  This  yields  curve  C  in  Fig. 
xx.  We  then  attempt  to  fit  elastic  jumps,  as  characterized  by 
equation  (18),  to  bridge  between  the  two  solution  curves  B  and  C. 
The  actual  location  of  the  elastic  jump  £  is  that  point  at  which 
the  jump  conditions  match  curves  B  and  C. 


The  location  of  the  FLS  and  the  maximum  flow  rates  for 
several  lung  volumes  are  shown  in  Fig.  5.  The  different  symbols 
correspond  to  several  different  combinations  of  friction  and 
stiffness  coefficients. 

Figures  6a  and  6b  illustrate  the  distribution  of  area  ratio 
(a)  and  speed  index  (S),  respectively,  along  the  continuous 
branching  network.  Each  curve  corresponds,  for  example,  to  a 
different  level  of  respiratory  effort;  those  curves  for  which  S 


either  reaches  or  exceeds  one  are  flow  limited.  Note  that  all 
area  ratios  are  less  than  one  (in  the  "collapsed"  portion  of  the 
tube  law) . 

The  critical  flow  curve  (corresponding  to  the  maximum  flow 
rate  for  given  lung  volume)  can  have  either  a  subcritical  or  a 
supercritical  branch  downstream  of  the  FLS ,  depending  on  pressure 
at  the  mouth  (£-1).  Supercritical  branches  either  remain 
supercritical  to  the  end  of  the  network  (£  -  1)  or  undergo  a 
transition  to  S<1  by  means  of  an  elastic  jump,  a  shock-like 
transition  from  supercritical  to  subcritical  flow  speed,  as 
described  above. 

Fig.  (6c)  shows  the  results  for  evacuation  of  a  lung  by 
applying  vacuum  at  the  mouth.  In  this  case  upstream  pressure  is 
held  fixed  and  mouth  pressure  is  reduced.  Pressure  progressivley 
falls  toward  the  mouth,  reaching  a  local  minimum  at  some  point 
downstream  of  the  FLS. 

The  pressure  curves  can  also  be  plotted  to  simulate  forced 
expiration  produced  by  respiratory  effort  by  simply  shifting  all 
the  curves  in  Fig.  (6c)  vertically,  so  that  P  ..-p  -0  as 

shown  in  Fig.  6(d)  (and  as  described  in  the  dSsHussI&S1? .  The 
latter  curves  represent  forced  expirations  performed  at 
different  levels  of  muscular  effort  producing  different  upstream 
conditions. 

An  important  feature  of  flows  with  elastic  Jumps  is  that  the 
amount  of  energy  dissipated  in  the  region  mouthward  from  the  FLS 
can  vary  depending  upon  the  location  of  the  jump.  This  occurs  in 
two  ways:  first,  the  frictional  dissipation  due  to  wall  shear 
stress  varies  depending  on  the  relative  lengths  of  supercrictical 
and  subcritical  flow;  and  second,  the  dissipation  associated  with 
the  jump  itself  depends  on  conditions  upstream  of  the  Jump. 

Thus,  it  becomes  possible  to  attain  different  pressures  at  the 
mouth,  for  identical  flows  upstream  of  the  FLS,  by  changing  the 
position  of  the  jump.  Reductions  in  mouth  pressure  cause  the 
elastic  Jump  to  move  progressively  mouthward  while  leaving  the 
flow  rate  exiting  from  the  lung  unchanged. 

Plotting  driving  pressure,  P  ,  -P  .  ,  from  the  different 
curves  of  Fig.  6(d)  against  the  corresponding  flow  rate,  at  fixed 
lung  volume,  yields  the  Iso-Volume  Pressure  Flow  (IVPF)  curves 
shown  in  Fig.  7.  The  critical  driving  pressure,  P.-P  ,  -P  .  , 
required  to  achieve  flow  limitation  is  indicated  on  e&cK  I^P? 
curve  in  Fig.  7  and  plotted  as  a  function  of  lung  volume  in  Fig. 

9  for  several  different  conditions. 

The  model  is  also  used  to  calculate  the  MEFV  curves  and  flow 
limitation  sites  for  different  gases  as  shown  in  Fig.  8. 


The  maximal  flow  rate  and  the  location  of  flow  limiting 
sites,  predicted  by  the  present  model,  are  summarized  in  Fig.  5, 
together  with  some  representative  experimental  data.  The 
agreement  is  generally  good  although  the  flow-volume  pattern 
obtained  by  simulation  appears  to  be  better  matched  to  the  canine 
experiments  than  the  human  experiments.  There  are  several 
possible  explanations  for  this.  First,  since  the  data  we  used  to 
produce  our  pressure-area  law  was  obtained  from  dogs,  we  might 
expect  the  predictions  to  exhibit  some  features  of  a  canine  flow- 
volume  curve.  Second,  if  the  difference  in  shape  of  the  MEFV 
curve  is  due  to  variations  in  tracheal  tension,  as  suggested  by 
Smaldone  and  Smith  (1985),  then  it  would  be  difficult  for  any 
model  not  incorporating  such  effects  to  distinguish  one  pattern 
from  the  other.  Yet  another  possibility,  related  to 
uncertainties  in  selecting  some  of  the  parameters  of  the  model, 
is  discussed  later. 

The  pressure  distributions  of  Fig.  6  are  of  particular 
interest  in  view  of  several  inferences  drawn  from  experimental 
measurements.  As  mentioned  in  the  Introduction,  it  is  generally 
held  that  the  FLS  is  the  furthest  mouthward  point  which  is 
insensitive  to  changes  in  downstream  pressure.  In  Fig.  6(c)  the 
FLS  is  identified  by  the  closed  circle  (£-0.56)  and  flow 
limitation  is  attained  for  mouth  pressures  more  than  59  cm  H„0 
below  alveolar  pressure.  As  mouth  pressure  falls  below  this2 
value,  the  region  within  which  pressure  is  insensitive  to  changes 
in  mouth  pressure  expands  some  distance  downstream  of  the  FLS, 
producing  a  region  of  supercritical  flow  (S>1)  terminated  by  an 
elastic  jump.  The  point  Identified  experimentally  as  the  FLS 
would,  according  to  these  results,  be  the  location  of  the  elastic 
Jump  rather  than  the  location  where  wave  speed  exceeds  flow  speed 
(S-l);lt  would  therefore  vary  depending  upon  the  values  of  mouth 
pressure  used  in  the  experiments. 

The  curves  in  Fig.  6  can  also  be  compared  to  the  area  and 
pressure  distributions  of  Lambert  and  Wilson  (LW)  (their  Fig. 
11).  Whereas  area  and  pressure  experience  a  relatively  smooth 
transition  through  the  FLS  in  our  calculations,  it  appears  that 
the  spatial  derivatives  of  both  area  and  pressure  become 
extremely  steep  (perhaps  approaching  infinity)  in  the 
calculations  of  LW.  This  is  at  least  suggestive  of  a  choke  point 
singularity  as  opposed  to  a  smooth  transition;  if  so,  it  would 
therefore  lead,  in  their  calculations,  to  a  solution  which  cannot 
be  mathematically  continued  downstream  of  the  FLS  and  which  is 
therefore  physically  impossible.  In  practical  terms,  it  would 
signify  that  they  had  overestimated  the  maximum  flow  possible. 

EcnilYalenoY  Between  Expiration  and  Evacuation 

During  a  natural  forced  expiration,  gas  is  driven  mouthward 
due  to  contraction  of  the  expiratory  muscles  and  the  muscles  of 
the  diaphragm  producing  pleural  pressures  much  in  excess  of 


atmospheric  pressure  At  each,  lung  volume,  alveolar  pressure 
exceeds  pleural  pressure  by  the  elastic  recoil  pressure  (pt?d)> 
which  depends  on  lung  volume .  In  experimental  investigations 
with  anesthetized  dogs  or  with  excised  lungs,  a  forced  expiration 
is  produced  by  evacuation  of  the  lung  by  means  of  a  vacuum 
applied  at  the  mouth  or  at  the  tracheal  exit.  That  these  two 
expiratory  maneuvers  are  essentially  equivalent,  can  be 
substantiated  by  the  following  line  of  reasoning. 

As  discussed  earlier,  the  flow  occuring  at  a  particular  lung 
volume  during  a  forced  expiration  is  equivalent  to  a  steady  flow 
at  the  same  lung  volume,  calculated  as  though  alveolar  volume 
remained  constant.  With  this  assumption  incorporated  into  the 
model,  the  flow  rate  through  the  bronchial  system  (and,  indeed, 
the  entire  solution)  is  found  to  depend  only  upon  the  transmural 
pressures  acting  at  the  upstream  (£-0)  and  downstream  (£-1)  end. 
Since  upstream  (peripheral)  transmural  pressure  is  fixed 
(P.  =P  -  at  £=0),  flow  rate  is  only  a  function  of  downstream 
traEsmural  pressure.  Noting  that  P  .  -P  ,• +P™  and 
(P,  ).»  ,-<dP  +P  .  ,  the  downstream  transmurax  pressure  can  be 
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pressure  drop  across  tSe  glottis  and  mouth;  a  const ant gvalue  for 
a  given  flow  and  lung  volume.  Since  fixed  values  for  this 
expression  yield  the  same  flow  rate,  it  can  be  seen  that  the 
effect  of  increasing  pleural  pressure  with  mouth  pressure  fixed, 
by  changing  expiratory  effort,  is  entirely  equivalent  to  the 
effect  of  holding  pleural  pressure  fixed  and  reducing  mouth 
pressure  below  atmospheric  by  means  of  vacuum  provided  the 
magnitude  of  pressure  change  in  each  instance  is  the  same. 

This  conclusion,  however,  is  not  necessarily  true  if  the 
critical  point  lies  in  the  region  where  the  "external"  pressure 
on  the  airway  varies  because  the  latter  passes  out  of  the 
thoracic  cavity. 

Breathing  Different  Gases 

Experimental  measurements  have  demonstrated  that  MEFV  curves 
differ  for  different  gases  at  high  and  moderate  lung  volumes,  but 
tend  to  merge  near  end  expiration  (Wood  et  al.  1976,  Mink  et  al., 
1979).  We  used  the  present  model  to  predict  flows  for  a  mixture 
of  80%  He  and  20%  02  and  a  mixture  of  80%  SFg  and  20%  0«.  The 
results  for  these  mixtures  and  air  are  shownDin  Fig.  8.  While 
the  general  trend  of  the  maximal  flow  curve  is  as  expected  (the 
difference  at  a  given  lung  volume  varies  as  the  square  root  of 
the  density) ,  a  somewhat  greater  density  dependence  is  predicted 
than  is  observed  experimentally,  particularly  at  low  lung 
volumes. 

These  differences  between  the  present  calculations  and 
experimental  resu3ts  may  be  due  to  a  stronger  dependence  on  gas 
viscosity  in  the  lung  than  predicted  by  the  model  and  may  be 
associated  with  the  differences  in  maximal  flow  at  low  lung 


volume  noted  earlier.  It  may  also  be  related  to  deficiencies  in 
the  tube  law  for  the  peripheral  airways,  where  there  is  a  paucity 
of  experimental  data. 

The  location  of  the  FLS  is  seen  also  to  depend  on  gas 
properties  and  tends  to  move  downstream  (mouthward)  when  gas 
density  is  reduced. 

The  Effect  of  Variations  in  External  Pressure 

Up  to  this  point  we  have  assumed  that  the  effective  external 
pressure  for  all  the  bronchial  airways  is  uniform  and  equal  to 
alveolar  pressure.  This  assumption  is  not  strictly  true  and  may 
in  fact  be  a  very  poor  approximation  in  the  extra-pulmonary 
bronchi  and  the  extra-thoracic  trachea.  It  is  possible,  however, 
to  relax  this  constraint  in  the  model  and  consider  the 
implications  of  variations  in  effective  external  pressure. 

To  examine  these  effects  we  change  the  model  and  assume  that 
all  lntrathoraclc  bronchi  are  exposed  to  the  same  uniform 
external  pressure  (extra-pulmonary  but  introthoracic  bronchi  are 
enclosed  by  the  lobes  of  the  lung  and  are  therefore  probably 
subjected  to  similar  loadings  as  the  intrapulmonary  airways).  We 
impose  a  drop  in  external  pressure  in  the  downstream  half  of  the 
trachea  as  shown  in  Fig.  10.  The  effective  plueral  pressure 
during  maximal  effort  is  taken  as  the  sum  of  maximal  dynamic 
effort  (Hyatt  and  Flath,  1966)  and  elastlo  recoil  pressure. 
However,  since  the  maximal  dynamic  effort  depends  on  the  flow 
rate  (the  result  of  the  calculation)  we  use  the  maximal  possible 
values  given  by  Hyatt  et  al.,  (1966),  as  representative  of  the 
greatest  potential  influence  of  an  external  pressure  gradient. 

The  loci  of  critical  points  for  maximal  lung  volume  (VL-100% 
TLC)  with  uniform  external  pressure,  and  with  a  presumed  gradient 
in  external  pressure  are  shown  in  Fig.  11.  The  differences  are 
seen  to  be  restricted  to  that  region  where  the  pressure  gradient 
is  non-zero.  In  Fig.  11  the  open  circles  signify  the  FLS 
locations  which  satisfy  the  upstream  boundary  conditions.  The 
results  for  maximal  flow  rate  and  the  predicted  flow  limiting 
sites  are  given  in  Fig.  12.  The  flow  rate  is  virtually 
unaffected  by  the  pressure  gradient  at  lower  lung  volumes  when 
the  FLS  is  located  further  upstream,  away  from  the  pressure 
gradient.  The  main  observable  effect  of  the  pressure  gradient  is 
to  displace  the  FLS  upstream  in  those  cases  for  which  the  FLS  is 
found  be  in  the  zone  of  a  significant  gradient  in  external 
pressure.  These  results  Justify  our  earlier  assumption  of 
uniform  external  pressure.  We  note,  however,  that  a  reduction  in 
external  pressure  acting  on  the  extrathoracic  airways  might 
significantly  alter  their  cross-sectional  area  even  though  the 
limiting  flow  rate  is  unchanged. 

As  noted  earlier,  the  gradient  of  effective  external 
pressure  may  negate  the  equivalence  between  a  forced  expiration 
and  a  forced  evacuation. 
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The  process  by  which  the  transition  point  is  reached  is 
described  in  the  section  on  the  model.  In  essence,  the  numerator 
of  the  RHS  of  equations  (16)  and  (17)  must  be  positive  for  S<1, 
pass  through  zero  when  S-l,  and  be  negative  when  S>1.  The 
magnitude  and  sign  of  the  numerator  depends,  in  general,  on  a 
complex  interaction  between  area  change  and  variations  in 
stiffness  and  friction.  Since  A  is  generally  decreasing  from 
£=0  to  $-1,  the  contribution  due°to  area  change  tends  to  drive 
the  flow  toward  S-l,  except  in  the  trachea  where  the  effect  of 
area  change  is  small.  According  to  earlier  assumptions,  vessel 
stiffness  monotonically  increases.  The  corresponding  term  in  the 
numerator  of  equation  (16)  and  (17),  as  S  approaches  1,  is 
therefore  negative,  driving  S  away  from  unity.  The  friction 
term,  which  always  retains  the  same  sign,  drives  S  toward  unity 
in  all  cases.  The  dominant  balance  in  this  equation  is  between 
friction  (driving  S  toward  1)  and  variations  in  stiffness 
(driving  S  away  from  1).  Consequently,  for  a  given  network 
geometry,  the  location  of  the  FLS  is  primarily  determined  by  the 
effects  of  friction  and  of  variations  in  stiffness  along  the 
conducting  system.  Increasing  the  stiffness  coefficient  (C.  )  or 
decreasing  the  friction  coefficient  (fT)  both  drive  the  FLS 
upstream.  Since  at  the  FLS  the  Reynolds  number  is  large 
(typically  greater  than  10*),  friction  is  dominated  by  the  second 
of  the  two  terms  in  equation  (15).  Therefore,  high  values  for  b 
drive  the  FLS  downstream;  and  in  fact,  for  b-0.0035  (the  value 
obtained  by  Reynolds  and  used  by  Wilson  and  Lambert  )  there  are 
no  transition  points  within  the  conducting  system  and.  Instead, 
choking  occurs  at  the  downstream  end  (£-1).  This  is 
characterlsitc  of  friction-dominated  flows  as  discussed  by 
Shapiro  (1977). 

For  the  values  used  in  this  study,  (b-0.0013;  C. -2. 40) the 
FLS  is  located  in  the  trachea  at  large  lung  volumes  and  moves 
into  the  second  generation  at  approximately  50%  of  TLC.  Under  no 
circumstances,  in  fact,  with  the  specified  parameters,  does  the 
FLS  move  further  upstream  than  the  second  generation.  With  other 
specified  distributions  of  rest  area,  tube  wall  properties  and 
frictional  behavior,  however,  this  result  too  can  vary.  In  the 
standard  simulation  the  laminar  friction  behavior  at  the 
peripheral  generation  is  such  that  the  frictional  dissipation  is 
1.5  (a  in  equation  (15))  times  as  large  as  that  in  Polseullle 
flow  (Reynolds,  1982).  The  excess  dissipation  is  presumably 
associated  with  the  effects  of  branching  and  curvature.  We 
found,  however,  that  the  results  of  the  simulation  are  virtually 
Independent  of  the  values  used  for  a;  in  doubling  a  to  3.0,  we 
obtained  results  which  are  indistinguishable  in  Fig  5. 

It  is  interesting  to  note  that  the  maximal  (critical)  flow 
rate  at  a  fixed  lung  volume  is  relatively  independent  of  either 
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the  frictional  term  or  changes  in  the  mechanical  properties  of 
the  conducting  system.  This  can  be  seen  in  Fig.  5a  in  which  all 
the  computed  flow  rates  for  a  system  with  fixed  geometry  but 
with  different  friction  and  stiffness  distributions  are  in  a 
narrow  band  with  deviations  of  less  than  10%. 

Physiological  Limitation 

The  experimental  data  and  the  predicitons  of  the  model 
presented  in  Figure  5(a)  exhibit  clear  differences  in  maximal 
flow  both  at  high  and  low  lung  volume.  The  differences  at  low 
lung  volume  have  been  discussed  above  in  the  context  of  the 
effects  of  gas  viscosity  and  will  be  addressed  again  later.  At 
high  lung  volume,  the  differences  are  caused  by  limitations 
imposed  by  the  respiratory  muscles  as  first  demonstrated  by  Hyatt 
et  al  .  In  this  case,  the  respiratory  muscles  are  Incapable  of 
producing  sufficient  force  to  reach  the  point  of  wavespeed  flow 
limitation.  Instead,  maximal  flow  is  observed  to  be  effort 
dependent . 

To  incorporate  the  effects  of  physiological  limitation  into 
the  model,  we  used  the  data  of  Hyatt  et  al  to  predict  the  maximal 
effort  associated  with  each  lung  volume.  Their  results  can  be 
plotted  as  in  Figure  7  in  terms  of  maximal  presssure  as  a 
function  of  the  rate  at  which  the  lung  is  emptying  (related  to 
the  velocity  of  muscle  contraction) .  The  intersection  of  these 
curves  with  the  maximal  flow  curve  predicted  by  the  model 
determines  the  true  maximal  flow  at  each  lung  volume.  Those  lung 
volumes  for  which  the  intersection  occurs  before  the  plateau  in 
the  flow  curve  is  reached  will  exhibit  effort  dependent  maximal 
flow.  Our  simulation  predicts  that  this  would  occur  at  a  lung 
volume  of  about  85%  TLC  as  compared  to  a  value  of  about  75% 
according  to  Mead  et  al.  Note  that  if  expiration  is  produced  by 
applying  progressively  lower  pressures  at  the  mouth,  then  the 
point  of  wavespeed  limitation  will  ultimately  be  attained, 
regardless  of  lung  volume  and  muscular  capability.  Accordingly, 
the  MEFV  curves  produced  by  forced  evacuation  will  differ  from 
those  fo  a  forced  expiration  at  high  lung  volumes,  a  difference 
which  has  been  observed  experimentally. 

Effects  of  a  Local  Perturbation 

In  our  previous  work  ,  we  investigated  the  effects  of 
applying  a  local  perturbation  in  the  distribution  of  airway 
stiffness,  reference  cross-sectional  area,  or  external  pressure. 
In  those  studies,  we  found  that  the  flow  can  be  driven  toward  a 
transition  point  by  any  of  the  following:  (i)  a  local  minimum  in 
Aq,  (il)  a  local  minimum  in  K_,  and  (iii)  a  local  maximum  in  the 
effective  external  pressure.  uIn  the  presence  of  one  of  these 
factors,  the  transition  to  supercritical  speed  might  occur  at  a 
more  peripheral  point  in  the  lung  when  either  normal  or 
pathological  variations  in  airway  structure  or  geometry  produce 
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any  of  the  above  conditions.  Moreover,  these  points  will 
continue  to  be  preferential  sites  for  flow  limitation  as  lung 
volume  decreases;  then,  rather  than  the  smooth  progress  of  the 
FLS  toward  the  periphery  predicted  in  the  present  calculations, 
the  FLS  might  progress  upstream  at  a  non-uniform  rate.  With  a 
sufficiently  strong  perturbation,  the  FLS  might  tend  to  jump  to 
or  from  tbe  region  of  the  perturbation  in  a  quasi-discontinuous 
fashion,  producing  an  observable  'glitch'  in  the  flow-volume 
curve . 


It  is  of  interest  to  note  that  the  model  developed  by  LW  was 
better  able  to  simulate  maximal  flows  and  the  effects  of  gas 
viscosity  at  low  lung  volumes  than  the  model  presented  here. 

They  accomplished  this  by  adjusting  their  pressure-area 
relationships  for  the  small  airways  until  the  predicted  flow 
matched  the  flows  obtained  experimentally  with  air.  In  doing  so, 
their  small  airways  were  made  extremely  compliant  (e.g.,  their 
16th  generation  airway  undergoes  a  ten-fold  increase  in  cross- 
sectional  area  as  transpulmonary  pressure  increases  from  zero  to 
3  cmH?0.  In  addition,  the  pressure-area  law  they  used  for 
negative  transmural  pressure  (encompassing  much  of  the  airway 
system  during  a  forced  expiration,  especially  at  low  lung 
volumes)  was  based  on  an  extrapolation  from  data  obtained 
exclusively  in  the  positive  pressure  range. 

In  the  LW  study,  as  in  ours,  the  greatest  weakness  lies  in 
the  lack  of  reliable  and  appropriate  experimental  data  on  which 
to  base  the  model.  While  the  tube  law  used  in  our  calculations 
is  more  closely  related  to  actual  experimental  data,  it  was 
necessary  to  make  several  assumptions  (most  notably  the  variation 
in  stiffness  in  the  airway  system)  that  are  only  weakly 
supported.  Another  potential  deficiency,  possibly  related  to  the 
failure  of  our  model  at  low  lung  volumes,  concerns  the  friction 
law  used  in  the  calculations.  This  is  essentially  the  same  as 
that  used  by  LW  but  with  different  values  for  the  two  constants. 
The  method  used  to  infer  a  local  pressure  drop  relationship  from 
measurements  across  an  entire  network  may  produce  a  misleading 
result  for  the  following  reason.  Consider,  for  example,  the 
behavior  at  the  extremes  of  high  and  low  Reynolds  number.  At 
sufficiently  low  Re,  the  flow  should  readily  adjust  at  each 
Junction.  Consequently,  the  velocity  profiles  will  be  nearly 
parabolic,  and  the  pressure  drop  will  be  very  nearly  that  of 
Pouseuille  flow  (giving  a-1.0).  At  high  values  of  Re,  energy 
dissipation  due  to  turbulence  will  dominate  and,  even  allowing 
for  the  loss  of  energy  due  to  the  generation  of  secondary  flow, 
the  pressure  drop  should  be  near  that  associated  with  a 
moderately  rough  pipe.  Both  estimates  are  considerably  below  the 
inferred  values  of  Reynolds  .  The  explanation  for  this  apparent 
discrepancy  may  be  associated  with  the  flow  behavior  at 
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intermediate  Renolds  numbers,  in  the  range  10 < Re <1000,  where  the 
first  term  in  equation  (15)  still  dominates.  Due  to  the  combined 
effects  of  the  Induced  secondary  motion  and  the  net  reduction  in 
flow  area  as  seen  by  gas  exiting  from  the  lung,  the  velocity 
profiles  will  be  quite  blunt  rather  than  parabolic.  Hence  the 
laminar  coefficient  in  equation  (15)  should  increase  with 
increasing  Re,  rather  than  remain  constant  as  it  is  forced  to  do 
in  the  current  formulation.  In  forcing  all  the  data  to  conform 
to  the  simple  form  of  equation  (15),  the  effects  of  an  increase 
in  viscous  shear  stress  (an  effective  increase  in  a)  is  reflected 
instead  by  an  increase  in  both  a  and  b.  Increasing  b  produces 
the  correct  overall  pressure  drop,  but  has  the  undesired  effect 
of  reducing  the  real  effect  of  gas  viscosity  as  predicted  by  the 
model . 

A  significant  increase  in  a  would  tend  to  Increase  the 
influence  of  gas  viscosity  on  maximal  flows  without  the  necessity 
of  malting  the  peripheral  airways  more  compliant.  However,  to 
determine  the  ultimate  validity  of  these  arguments  will  require 
more  definitive  data  both  on  the  local  pressure  drop 
characteristics,  and  on  the  constitutive  relations  for  the 
smaller  airways. 


The  model  presented  here  represents  one  step  toward 
developing  a  more  realistic  simulation  of  a  forced  expiration. 

We  want  to  emphasize,  though,  the  need  for  further 
investigations.  In  this  initial  study,  we  have  identified 
several  areas  where  additional  research  is  needed  before  the 
models  can  be  refined.  These  include:  (1)  the  tube  law, 
especially  for  the  peripheral  airways,  (2)  the  expression  used  to 
describe  the  pressure  drop  due  to  frictional  effects,  (3)  the 
behavior  of  supercritical  flow  and  elastic  jumps  in  pulmonary 
airways,  and  (4),  the  effects  of  asymmetry  in  the  branching 
geometry  and  distribution  of  flow.  The  first  two  have  been 
discussed  previously.  The  third  point  arises  from  the 
recognition  that  nearly  all  studies  of  supercritical  flow  and 
elastic  Jumps  have  been  conducted  in  latex  rubber  (penrose)  tubes 
filled  with  water.  It  is  entirely  possible  that  the  stability  of 
such  flows  is  different  from  that  in  the  flow  of  gas  through  the 
conducting  airways.  The  presence  of  instabilities  could  lead  to 
solutions  of  a  very  different  character  downstream  of  the  FLS 
from  those  predicted  by  the  model  since  the  latter  presumes  the 
existence  and  stability  of  a  quasi-steady  flow.  And  finally,  it 
is  essential  that  we  recognize  the  potentially  Important 
implications  of  asymmetry  .  While  it  is  gratifying  to  see  the 
degree  of  agreement  between  this  simple  model  and  real  flows,  it 
is,  at  some  time,  surprising  that  a  structure  so  complex  yields 
to  so  simple  an  analysis.  There  is  clearly  much  to  be  learned 
from  the  effects  of  lung  asymmetry. 


IV.  PATIENT  STUDIES 


One  of  the  outcomes  of  our  model] ing  efforts  was  a  better  appreciation 
of  how  maximum  expiratory  flow  could  influence  lung  gas  volume  during  HFV.  In 
our  progress  report  last  year,  we  described  experiments  in  which  the  mechanism 
of  the  hyperinflation  observed  during  HFV  was  thought  to  be  expiratory  flow 
limitation.  In  this  year  we  describe  additional  experiments  designed  to 
approach  this  hypothesis.  In  order  to  make  appreciation  of  both  sets  of 
experiments  simpler,  we  will  include  both  sets  of  results  in  this  year's 
progress  report. 
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1 .  Introduction. 

It  has  now  been  well  established  that  dynamic  hyperinflation  of  the 
lungs  may  occur  in  animals  (Simon  et  al,  1984)  or  in  patients  ( Reheder  and 
Didier,  1984;  Saari  et  al,  1984)  during  low  tidal  volume,  high  frequency 
oscillatory  ventilation  (HFOV).  Simon  et  al  (1984)  originally  attributed  this 
phenomenon,  which  they  observed  in  normal  dogs,  to  asymmetry  of  pulmonary 
impedance  during  inspiration  and  expiration.  We  have  noted,  based  on  analysis 
of  cineradiobronchographs  ( Gavriely  et  al,  1985),  that  asymmetric  variations  in 
airway  diameter  about  the  mean  occur  in  canine  central  airways  during  HFOV. 

When  this  occurs,  airway  opening  pressure  tracings  show  similar  asymmetry,  with 
larger  negative  deflections  than  positive  deflections  from  the  mean.  Similar 
airway  opening  pressure  asymmetries  have  been  observed  in  patients  during  HFOV. 

These  studies  were  undertaken  in  patients  with  respiratory  failure  to 
address  the  mechanisms  responsible  for  dynamic  hyperinflation  and  the  associ¬ 
ated  airway  opening  pressure  fluctuation  asymmetry.  To  do  so,  we  measured 
changes  in  lung  volume  during  HFOV,  and  compared  these  results  with  steady 
maximum  expiratory  flow  curves,  as  well  as  with  predictions  made  by  a  series  of 
computational  models,  each  of  which  incorporated  a  different  mechanism  which 
could  lead  to  dynamic  hyperinflation.  We  found  that  expiratory  flow  limit¬ 
ation,  rather  than  fixed  asymmetric  or  variable  asymmetric  resistance,  best 
accounted  for  our  experimental  observations. 


2 .  Methods. 

Patient  Studies 

Five  female  patients  (ages  28-86)  with  chronic  respiratory  insufficiency 
served  as  subjects  for  the  study.  Each  required  long-term  mechanical  ventila¬ 
tory  support  for  underlying  neuromuscular  disease  (patients  #1-4)  or  for  severe 
obstructive  lung  disease  (patient  #5).  Informed  consent  for  the  protocol  as 
institutionally  approved  was  obtained  from  each  patient.  The  patients 
tolerated  applications  of  HFOV  (see  below)  well,  without  apparent  or  stated 
discomfort. 

Changes  in  lung  volume  were  assessed  with  a  Respitrace  impedance 
plethysmograph.  Belts  were  positioned  around  the  abdomen  and  around  the  chest; 
relative  gains  were  adjusted  so  that  changes  in  the  sum  of  the  two  signals  were 
directly  proportional  to  changes  in  mean  lung  volume.  Calibration  of  the 
Respitrace  volume  signal  was  performed  by  delivering  a  number  of  increments  of 
gas  from  a  calibrated  syringe  into  the  patient.  The  measured  volume  signal  was 
linear  over  the  range  of  volumes  studied;  the  calibration  was  unaffected  by 
applied  oscillation  as  demonstrated  previously  (Saari  et  al,  1984).  Each 
patient  was  studied  in  the  supine  position  and  did  not  change  posture  or  body 
position  during  the  study. 

Two  series  of  experiments  were  performed.  In  Series  A,  the  dependence 
of  lung  volume  on  HFOV  frequency  and  mean  airway  opening  pressure  was  deter¬ 
mined  in  three  subjects.  In  Series  B,  we  compared  maximal  expiratory  flow 
curves  with  the  peak  expiratory  flow  vs.  lung  volume  relationship  found  during 
HFOV  in  four  subjects. 


Series  A 

The  high  frequency  oscillatory  ventilator  included  a  sealed  piston 
oscillator  and  a  high  impedance  bias  flow  ( Rossing  et  al,  1984).  Air  or  an 
air-oxygen  mixture  was  bled  into  the  system  from  a  high  pressure  source  across 
a  needle  valve.  The  bias  outflow  (rate  16-20  L/min,  measured  with  a  calibrated 
rotameter)  was  servo-controlled  to  maintain  mean  airway  opening  pressure  at  a 
predetermined  level  (2-13  cm  H2O)  with  a  3  second  time  constant,  independent 
of  whether  or  not  oscillations  were  applied.  Airway  opening  pressure  (Pao)  was 
measured  through  a  lateral  side  arm  at  right  angles  to  the  13  nun  i.d.  HFOV 
delivery  tube,  using  a  solid  state  pressure  transducer  (Ailtech  MS-10);  mean 
airway  opening  pressure  was  obtained  electronically.  Fifty  or  100  ml 
sinusoidal  oscillations  (frequency  0.5-20  Hz)  were  provided  by  the  sealed 
piston  pump  driven  by  a  linear  magnetic  motor.  Tidal  volume  was  measured  by 
integrating  airway  flow  (measured  with  a  number  2  Fleisch  pneumotachograph  and 
a  Validyne  MP-45  +  50  cm  H2O  pressure  transducer)  placed  between  the  bias 
flow  inlet  and  the  airway  opening.  To  calibrate  the  tidal  volume  signal,  we 
used  the  HFOV  apparatus  to  deliver  tidal  volumes  into  and  out  of  a  large  rigid 
vessel  of  known  adiabatic  compliance.  Integrated  flow  was  related  to  delivered 
tidal  volume  (derived  from  pressure  changes  in  the  rigid  vessel)  at  each 
frequency-tidal  volume  combination  to  be  employed. 

The  patient's  tracheostomy  tube  was  disconnected  from  the  conventional 
ventilator  and  the  patient  allowed  to  exhale  to  FRC  ( Pao=0  cm  H2O)  for  5-10 
seconds,  until  the  patient's  lung  volume,  as  measured  by  the  Respitrace,  was 
stable.  The  subject  was  then  connected  to  the  HFOV  apparatus  with  the 
oscillator  turned  off,  at  which  time  the  lungs  were  inflated  by  the  bias  flow 
circuit  until  a  new  stable  lung  volume,  determined  by  the  prechosen  mean  Pao, 
was  reached.  Sinusoidal  oscillations  were  applied  at  a  prechosen  frequency  and 
tidal  volume  for  30-45  seconds;  further  changes  in  lung  volume  were  monitored 
by  the  impedance  plethysmograph,  and  the  subject  was  then  returned  to  the 
conventional  ventilator.  The  entire  sequence  was  repeated  as  HFOV  frequency 
and  tidal  volume  and  mean  airway  opening  pressure  were  systematically  varied. 

Series  B 

Maximal  expiratory  flow-volume  curves  were  generated  as  follows.  Each 
subject  was  allowed  to  relax  to  FRC;  the  patient's  tracheostomy  cannula  was 
then  connected  to  a  3  liter  syringe.  1500  ml  were  delivered  from  the  syringe 
into  the  patient;  2  seconds  later,  the  syringe  plunger  was  rapidly  and  forcibly 
withdrawn,  effecting  a  maximal  expiration.  Two  to  five  such  maximal  expira¬ 
tions  were  performed  over  a  10-15  minute  period.  Lung  volume  versus  time 
curves  during  these  maneuvers  (obtained  from  impedance  plethysmograph  tracings) 
were  averaged  as  volume  as  a  function  of  time,  and  a  maximum  expiratory 
flow-volume  curve  was  obtained  by  discrete  differentiation  of  the  average 
volume-time  curve  every  0.1  sec. 

After  completion  of  the  above  maneuvers,  the  HFOV  piston  pump  used  in 
Series  A  was  used  to  deliver  sinusoidal  oscillations  of  volume  50  to  400  ml  at 
frequencies  between  0.33  and  10  Hz.  The  bias  flow  apparatus  was  replaced  by  an 
oxygen-filled  weighted  spirometer,  which  supplied  gas  through  a  high  impedance 
to  oscillatory  flow  ( Saari  et  al,  1984)  such  that  mean  airway  opening  pressure 
was  maintained  at  2-3.5  cm  H2O.  Periods  of  HFOV  were  applied  as  in  Series  A 
until  a  new  stable  lung  volume  was  reached,  several  HFOV  frequency-tidal 
volume  combinations  were  applied  in  sequence  without  return  to  FRC  or 


mechanical  ventilation,  but  never  for  more  than  60  seconds  without  returning  to 
conventional  ventilation-  Changes  in  lung  volume  were  recorded  using  the 
impedance  plethysmograph. 


Computational  Model 

Three  electrical  circuits  (Figure  13)  were  used  to  model  the  mechanical 
behavior  of  the  high  frequency  ventilator  and  the  respiratory  system.  In  the 
analogy,  current  was  used  to  model  volume  flow,  and  voltage  was  used  to  reflect 
pressure  (atmospheric  pressure  was  considered  as  ground  potential).  The  HFOV 
pump  was  thus  a  current  source  (Ipump);  the  bias  flow  was  modelled  with  a  volt¬ 
age  source  (Vmap),  whose  value  represented  mean  airway  pressure,  coupled  to  the 
airway  opening  node  by  a  large  inductance  (Lbf).  This  configuration  resulted 
in  a  high  impedance  to  oscillatory  flow  (Butler  et  al,  1980).  The  respiratory 
tree  was  modelled  as  a  central  airway  "shunt"  compliance  (Caw),  which  reflects 
potential  radial  oscillatory  motion  of  central  airway  walls,  and  the  remaining 
lung  as  a  path  from  the  airway  opening  node  to  the  alveolar  volume,  modelled  as 
a  pure  compliance  (Calv).  Inspiratory  and  expiratory  limbs  of  the  latter  path 
were  isolated  by  diodes  to  facilitate  implementation  of  asymmetry  (Simon  et  al, 
1 984) . 


The  inspiratory  limb  in  all  three  models  consisted  of  a  fixed  inspira¬ 
tory  resistance  (Rinsp)  leading  to  the  alveolar  node.  The  structure  of  the 
expiratory  limb  was  varied  to  model  each  of  three  types  of  asymmetric  airways 
impedance.  Fixed  asymmetric  resistance  (Figure  la)  was  modelled  using  a  linear 
resistor  (Rexp);  variable  asymmetric  resistance  (Figure  1b)  was  modelled  with  a 
variable  resistance  whose  value  varied  inversely  with  lung  volume  [Rexp( Palv) ] ; 
and  expiratory  flow  limitation  (Figure  1c)  was  modelled  using  a  PNP  transistor 
to  limit  flow  (current)  in  series  with  a  resistance  (Rexp)  whose  value  was 
equal  to  that  of  Rinsp.  The  transistor  collector-emitter  voltage-current 
characteristic  (Gray  et  al,  1969)  was  chosen  to  mimic  the  isovolume 
pressure-flow  relationships  of  human  lungs  (Hyatt  et  al,  1958).  The  maximum 
allowable  expiratory  flow  through  the  transistor  is  determined  by  Palv  and 
Vcrit.  Vcrit  represent  the  value  of  Palv  at  which  there  is  no  allowable 
expiratory  flow.  As  Palv  (and  lung  volume)  rises  above  Vcrit,  the  maximum 
allowable  expiratory  flow  rises  proportionally  to  the  difference  between  Palv 
and  Vcrit.  Component  values  were  chosen  to  reflect  physiologic  parameters 

typical  of  patients  such  as  those  in  our  study,  and  are  listed  for  each  model 

in  Table  1.  Tidal  volume  was  set  to  100  ml;  frequencies  of  1  to  8  Hz  were 

used.  Mean  airway  pressures  of  3  and  10  cm  H2O  were  studied  for  each  type  of 

asymmetry. 

Simulations  of  HFOV  by  these  models  were  performed  by  an  iterative 
integration  computer  program  (SPICE2,  copyright  Univ.  of  Calif.,  Berkeley)  for 
a  sufficiently  long  interval  to  attain  a  steady  state.  Excursion  extremes  of 
Palv  and  Pao  were  taken  from  the  last  second  of  simulation,  and  mean  Palv  was 
estimated  as  the  average  of  excursion  extremes. 

3.  Results. 

Patient  Studies 
Series  A 

The  dependence  of  dynamic  hyperinflation  on  frequency  and  on  mean  airway 
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opening  pressure  is  illustrated  in  Figure  14.  When  mean  airway  opening 
pressure  was  low  (2-5  cm  H2O),  dynamic  hyperinflation  occurred  during  HFOV  at 
relatively  low  frequencies,  and  lung  volume  increased  with  frequency  such  that 
over  a  one  liter  dynamic  change  in  FRC  was  observed  in  2  subjects.  In 
contrast,  when  the  mean  airway  pressure  was  raised  (7-13  cm  H2O),  a  different 
dependence  of  dynamic  hyperinflation  on  frequency  was  observed  at  the  same 
Vf.  The  elevated  mean  airway  opening  pressure  was  associated  with  an 
elevated  mean  lung  volume  over  resting  FRC  before  oscillations  were  applied, 
and  there  was  no  further  increase  in  lung  volume  as  frequency  was  raised,  until 
the  frequency  reached  some  critical  value  (approximately  2-4  Hz).  Above  this 
critical  frequency,  lung  volume  increased  with  frequency  in  a  fashion  nearly 
parallel  to  that  observed  at  the  lower  mean  airway  opening  pressure^.  As 
demonstrated  previously,  the  degree  of  dynamic  hyperinflation  was  greater  at 
each  frequency  when  larger  tidal  volumes  were  applied  at  constant  mean  airway 
opening  pressure  (Figure  15). 

To  define  the  importance  of  expiratory  flow  rate  in  dynamic 
hyperinflation  during  HFOw,  we  varied  the  inspiratory-expiratory  duration  (I:E) 
ratio  of  the  oscillation  from  1:2  to  2:1  in  two  subjects,  during  HFOV  with 
fixed  V-j-  (100  ml)  and  mean  Pao  (2  cm  H2O).  A  quasi-sinusoidal  volume 
waveform  was  used.  Dynamic  hyperinflation  was  most  pronounced  at  any  frequency 
when  the  greatest  I:E  ratio  was  applied.  These  data  are  expressed  as  a 
function  of  mean  expiratory  flow  in  Figure  16. 

Series  B 

Figure  17  presents  expiratory  flow  versus  lung  volume  relationships 
found  during  steady  forced  maximal  expiratory  flow  and  during  HFOV 
applications.  HFOV  data  are  expressed  as  peak  expiratory  flow  at  the  airway 
opening  versus  lung  volume.  The  peak  expiratory  flow  during  HFOV  was  almost 
always  greater  than  that  found  during  the  steady  forced  maximal  expiration 
maneuver  at  the  same  lung  volume.  Subject  #1 ,  who  demonstrated  the  least 
difference  between  oscillatory  and  steady  expiratory  flow-volume  relationships, 
demonstrated  the  least  impaired  steady  maximal  expiratory  flow-volume  curve. 

In  addition,  she  was  the  youngest  patient,  and  clinically  had  the  least  severe 
functional  pulmonary  Impairment,  as  judged  by  usual  arterial  blood  gases, 
frequency  of  respiratory  infections,  and  chronic  inspired  oxygen  fraction 
requirement. 

Model  Predictions 

The  dependence  of  lung  volume  on  frequency  and  mean  airway  opening  pres¬ 
sure  predicted  by  each  of  the  three  asymmetric  respiratory  impedance  models  is 
illustrated  in  Figure  18.  In  the  first  model,  the  expiratory  limb  of  the 
respiratory  path  was  composed  of  a  fixed  resistance  whose  magnitude  was  twice 
that  of  the  inspiratory  resistance.  At  low  mean  Pao  (3  cm  H2O),  lung  volume 
progressively  increased  with  frequency  throughout  the  range  of  frequencies 
simulated.  Raising  mean  Pao  to  10  cm  H2O  raised  lung  volume  but  had  no 
effect  on  the  shape  of  the  lung  volume  -  frequency  relationship.  There  was  no 
range  of  frequencies  over  which  lung  volume  was  independent  of  frequency  at 
either  mean  Pao.  The  second  model  incorporated  an  expiratory  resistance  whose 
magnitude  varied  inversely  with  lung  volume.  At  low  mean  Pao,  there  was  a  more 
rapid  rise  of  lung  volume  with  frequency  than  was  seen  in  the  fixed  asymmetry 
model,  reflecting  the  greater  degree  of  impedance  asymmetry  at  low  lung  volumes 
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(Table  1);  however,  the  dependence  of  lung  volume  change  on  frequency  was 
otherwise  similar  to  that  found  for  the  case  of  fixed  asymmetry.  At  high  mean 
Pao,  the  variable  Rexp  model  predicted  somewhat  less  dynamic  hyperinflation  at 
the  higher  frequencies  than  did  the  fixed  asymmetry  model,  reflecting  the 
smaller  degree  of  asymmetry  at  high  lung  volumes  (Table  1),  but  there  was  still 
a  progressive  rise  of  lung  volume  with  frequency.  Thus,  predictions  of  both 
the  fixed  and  variable  expiratory  resistance  models  are  qualitatively  different 
from  the  dependence  of  lung  volume  on  HFOV  frequency  and  mean  Pao  observed  in 
our  Series  A  patients,  in  that  neither  demonstrates  a  plateau  of  lung  volume 
over  a  range  of  low  frequencies. 

In  contrast,  the  dependence  of  lung  volume  on  frequency  and  mean  Pao 
predicted  in  the  expiratory  flow  limitation  model  did  qualitatively  resemble 
that  observed  in  our  patients.  A  monotonic  increase  in  lung  volume  with 
frequency  was  evident  at  low  mean  Pao;  when  mean  Pao  was  raised,  a  plateau  of 
lung  volume  over  a  range  of  lower  frequencies  appeared. 

While  the  dependence  of  the  lung  volume  -  frequency  relationship  on  mean 
Pao  varied  substantially  among  the  three  models,  several  features  of  the 
models'  predictions  were  strikingly  similar.  First,  all  three  models  predicted 
airway  opening  pressure  fluctuations  with  qualitatively  similar  asymmetries 
about  their  means;  further,  the  diminution  of  fluctuation  asymmetry  with  time 
recently  reported  by  Rehder  and  Didier  (1984)  occurred  with  all  three  modes  of 
impedance  asymmetry.  Second,  the  time  delay  between  excursions  of  airway 
opening  (Pao)  and  alveolar  pressures  (Palv)  in  each  model  is  about  one-quarter 
cycle,  with  Pao  leading  Palv  (Figure  19).  Since  Pao  and  Palv  are  each 
proportional  to  the  tidal  volumes  delivered  to  the  central  airway  shunt 
compliance  and  alveolar  zones,  respectively,  these  two  parts  of  the  tidal 
volume  introduced  into  the  trachea  must  also  be  out  of  phase.  Thus,  the  sum  of 
the  "shunted"  and  "peripheral"  portions  of  tracheal  tidal  volume  is  greater 
than  the  tidal  volume  introduced  into  the  trachea  itself.  Figures  20a  and  20b 
show  that  the  magnitude  of  shunted  and  peripheral  tidal  volumes  are  closely 
similar  for  each  of  the  three  models  of  asymmetric  lung  impedance;  substantial 
"shunt"  airway  tidal  volumes  are  evident  at  the  higher  frequencies  in  each 
case.  Dependence  of  oscillatory  flow  delivered  to  the  alveoli  on  frequency  and 
mean  Pao,  shown  in  Figure  21,  is  nearly  identical  among  the  three  models  of 
lung  asymmetry.  Distal  oscillatory  flow  increased  roughly  linearly  with 
frequency  at  low  frequencies,  then  rose  less  rapidly  with  higher  frequencies 
above  a  critical  frequency,  about  4  Hz.  In  addition,  oscillatory  flow  in  the 
alveolar  zone  depended  little  on  mean  Pao. 

4.  Discussion. 

Dynamic  pulmonary  hyperinflation  has  been  shown  to  occur  during  high 
frequency,  low  tidal  volume  oscillatory  ventilation  (HFOV)  in  both  animals  and 
man.  Various  investigators  have  suggested  that  this  phenomenon  stems  from  dif¬ 
fering  inspiratory  and  expiratory  impedances  to  airflow.  Such  differences 
could  arise  from  the  asymmetries  in  resistive  pressure  dissipation  that  are 
known  to  occur  between  inspiration  and  expiration  in  a  rigid  branching  network 
(Reynolds  and  Lee,  1981).  Lung  volume  dependence  of  resistance  (Mead  and 
Whittenberger,  1956)  could  also  contribute  to  dynamic  hyperinflation,  and 
dynamic  narrowing  of  the  airways  reflecting  expiratory  flow  limitation  is  a 
third  possible  mechanism  contributing  to  dynamic  hyperinflation.  We  thought 


that  it  might  be  possible  to  distinguish  among  these  mechanisms  by  modelling 
the  electrical  analog  of  each  and  comparing  general  aspects  of  predictions  made 
from  each  model  to  data  sets  obtained  in  patients  demonstrating  hyperinflation 
during  HFOV.  Our  data  are  compatible  with  the  hypothesis  that  hyperinflation 
arises  as  a  result  of  expiratory  flow  limitation. 

Our  models  utilized  both  well  established  and  novel  electrical  analogs 
to  simulate  various  aspects  of  the  mechanical  behavior  of  the  respiratory  sys¬ 
tem.  All  three  models  incorporated  a  shunt  capacitor  (Caw)  to  model  central 
airway  compliance,  as  suggested  by  Head  (1969).  During  an  oscillatory  cycle, 
this  capacitor  receives  a  fraction  of  the  tidal  volume  delivered  at  the  airway 
opening,  and  this  fraction  varies  with  the  oscillatory  frequency.  Although  up 
to  70f  of  the  airway  opening  tidal  volume  was  delivered  onto  the  airway  shunt 
capacitor  (at  high  frequency),  this  volume  was  then  exhaled  into  the  alveolar 
zone  compliance.  Thus,  to  the  extent  that  our  models  predict  real  lung 
behavior,  central  airway  shunting  might  decrease  HFOV  efficiency  somewhat  by 
absorbing  a  fraction  of  the  tidal  volume  delivered  at  the  airway  opening 
( Rossing  et  al. ,  1981).  However,  this  mechanism  might  also  serve  to  enhance 
mixing  between  the  central  airways  and  the  alveolar  zone. 

A  number  of  novel  features  were  incorporated  into  our  electrical 
analogs.  First,  we  isolated  the  inspiratory  and  expiratory  limbs  of  the  model 
using  idealized  diodes,  as  suggested  recently  by  Simon  and  coworkers  (1984). 
This  allowed  independent  manipulation  of  the  expiratory  limb  without  alteration 
of  the  inspiratory  limb  or  the  rest  of  the  model.  Second,  we  simulated  the 
role  of  the  bias  flow  in  maintaining  constant  mean  airway  opening  pressure  by 
using  a  voltage  (pressure)  source  connected  to  the  airway  opening  by  a  large 
inductance.  Since  the  impedance  of  this  inductor  increases  with  frequency,  it 
excludes  oscillatory  flow  while  allowing  steady  mean  airway  opening  pressure  to 
be  maintained,  after  an  initial  transient.  Note  that,  as  during  HFOV,  the  bias 
flow  pressure  source  provides  the  volume  which  enters  the  lung  during  hyperin¬ 
flation;  it  is  the  only  source  of  volume  (charge)  in  these  models.  Third,  we 
have  employed  a  bipolar  transistor  as  a  lumped  parameter  model  of  expiratory 
flow  limitation.  In  the  configuration  used,  the  transistor  allows  flow  through 
its  axial  path  of  magnitude  less  than  or  equal  to  a  maximum  value.  This 
maximum  flow  limit  is  proportional  to  the  pressure  difference  between  Palv  and 
Vcrit  (the  critical  airway  closing  pressure),  and  is  independent  of  the 
downstream  pressure,  provided  that  the  downstream  pressure  is  sufficiently  low 
to  cause  flow  limitation  to  occur.  The  only  way  to  increase  maximum  flow  in 
the  expiratory  limb  is  for  Palv  to  increase. 

Similar  mechanisms  also  account  for  dynamic  hyperinflation  in  the  other 
two  models;  differences  stem  only  from  the  determinants  of  the  degree  of 
asymmetry  in  each  model.  In  the  model  with  lung  volume  dependence  of  expira¬ 
tory  resistance,  the  relationships  between  Palv  and  frequency  for  the  two 
values  of  mean  airway  opening  pressure  examined  were  not  parallel;  this  results 
because  Rexp  is  lower  at  all  values  of  frequency  for  the  higher  mean  Pao,  and 
the  incremental  decrease  in  Rexp  for  unit  changes  in  lung  volume  are  less  at 
higher  lung  volumes.  Clearly,  the  precise  relationships  obtained  between  Palv 
(lung  volume)  and  frequency  at  various  values  of  mean  Pao  in  this  model  will 
depend  on  the  relationship  chosen  for  Rexp  as  a  function  of  Palv.  We  chose  an 
inverse  relationship  for  this  model  based  on  the  data  Dubois  and  co-workers 
(1956)  obtained  during  panting,  Of  course,  during  external  forced  oscillation 
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at  these  tidal  volumes  and  frequencies,  other  relationships  may  apply. 
Nevertheless,  both  models  which  do  not  incorporate  dynamic  flow  limitation 
predict  a  dependence  of  Palv  on  frequency  at  various  mean  airway  pressures 
substantially  different  from  that  which  we  observed  experimentally.  Thus,  it 
seems  likely  that  these  two  models  fail  to  incorporate  vital  features  of  the 
mechanical  behavior  of  our  patients’  respiratory  systems. 

In  contrast,  our  model  which  incorporated  dynamic  flow  limitation 
demonstrated  a  dependence  of  Palv  on  frequency  which  was  quite  similar  to  that 
observed  in  all  three  of  the  subjects  we  studied.  The  critical  feature  which 
was  uniquely  similar  in  both  the  data  and  the  flow  limitation  model  was  the 
relationship  between  Palv  and  frequency  at  low  frequencies  (less  than  2-3  Hz) 
at  high  and  low  mean  airway  opening  pressures.  This  results  from  the  fact  that 
the  maximum  expiratory  flow  allowed  through  the  transistor  depends  only  on  lung 
volume.  With  high  mean  Pao,  the  elevated  lung  volume  allows  a  maximum 
expiratory  flow  which  is  greater  than  that  required  at  low  frequencies.  In 
this  circumstance,  there  is  no  flow  limitation,  and  the  transistor  appears  as  a 
short-circuit,  with  no  resistance.  Thus,  at  frequencies  below  the  critical 
frequency  (3  Hz  for  mean  Pao  =  10  cm  H2O  and  Vf  =  100  ml),  there  is  no 
symmetry  of  inspiratory  and  expiratory  paths,  and  no  dynamic  hyperinflation 
occurs.  In  contrast,  with  lower  mean  Pao  (and  lower  initial  lung  volume),  the 
expiratory  flow  required  even  at  lower  frequencies  may  exceed  that  allowed  at 
the  lower  lung  volume. 

In  this  flow  limitation  model,  we  used  transistor  characteristics  which 
reflect  features  of  isovolume  pressure-flow  curves  obtained  during  steady 
flow.  We  realize  that  these  relationships  may  not  prevail  during  oscillatory 
flow,  and  therefore  compared  the  expiratory  flow-volume  relationships  during 
HFOV  to  those  during  forced  maximal  exhalations  in  Series  B.  Our  results 
clearly  demonstrate  that  oscillatory  and  steady  expiratory  flow-volume  relation¬ 
ships  need  not  be  the  same,  a  result  considered  in  more  detail  below.  Neverthe¬ 
less,  the  similarity  of  the  model’s  predictions  and  our  experimental  observa¬ 
tions  provides  support  for  the  hypothesis  that  expiratory  flow  limitation 
accounts  in  large  part  for  the  dynamic  hyperinflation  observed  during  HFOV. 

Further  support  for  the  hypothesis  that  dynamic  hyperinflation  results 
from  expiratory  flow  limitation  can  be  gathered  from  Figure  16.  By  varying 
inspiration-expiration  duration  ratio  (I:E),  we  varied  mean  expiratory  flow 
while  maintaining  tidal  volume,  frequency,  and  mean  airway  opening  pressure 
constant.  Though  I:E  ratio  substantially  influenced  lung  volume  change  for 
fixed  frequency,  tidal  volume,  and  mean  airway  opening  pressure,  its  had  little 
influence  when  change  in  lung  volume  was  considered  as  a  function  of  mean 
expiratory  flow.  This  result  indicates  that  changes  in  I:E  ratio  influenced 
the  degree  of  dynamic  hyperinflation  simply  by  changing  expiratory  flow,  which 
likely  determines  the  degree  of  dynamic  hyperinflation. 

In  Series  B,  we  directly  compared  the  lung  volumes  needed  to  sustain 
expiratory  flows  during  HFOV  with  the  flow-volume  relationship  observed  during 
a  steady  forced  maximal  expiration  (MEFV)  from  1.5  L  above  FRC.  At  any  lung 
volume,  peak  expiratory  flow  during  HFOV  exceeded  that  during  the  MEFV  maneuver 
by  several-fold.  A  number  of  phenomena  can  result  in  expiratory  flow  at  a 
given  lung  volume  greater  than  that  observed  during  a  steady  forced  maximal 
expiration.  First,  dynamic  collapse  of  central  airways  can  lead  to 


augmentation  of  expiratory  flow  at  the  onset  of  forced  expiration  (Knudson  3t 
al. ,  1974).  This  occurs  only  once  during  a  MEFV  maneuver  but,  during  HFOV, 
could  occur  with  every  cycle.  Such  repeated  airway  collapse  would  result  in 
far  greater  peak  expiratory  flow  at  the  trachea  than  out  of  the  alveoli;  the 
MEFV  curve,  though,  is  vastly  more  dependent  on  the  latter.  Second,  the 
distribution  of  oscillatory  flow  may  be  different,  at  any  lung  volume,  from 
that  of  maximal  expiratory  flow.  For  example,  during  a  partial  forced 
expiratory  flow  maneuver  (Melissinos  et  al. ,  1979),  initial  flow  favors  lung 
regions  with  short  time  constants;  it  may  be  that  flow  during  HFOV  originates 
from  the  "faster"  expiratory  pathways  of  a  non-homogeneous  lung.  Indeed,  in 
patients  with  obstructive  lung  disease,  the  expiratory  flow  during  a  partial 
forced  expiration  may  exceed  that  during  a  MEFV  maneuver  from  TLC  by  over  2 
times  (Takishima,  et  al. ,  1967).  Note  that  in  each  of  the  cases  above,  the 
mechanism  leading  to  dynamic  hyperinflation  is  still  expiratory  flow 
limitation,  even  though  a  different  expiratory  flow  versus  lung  volume 
characteristic  may  determine  the  degree  of  dynamic  hyperinflation. 

In  summary,  through  the  combined  use  of  electrical  analog  modelling  and 
experimental  observations  in  intubated  ventilator-dependent  patients,  we  pro¬ 
pose  that  expiratory  flow  limitation  is  responsible  for  the  impedance  asymmetry 
that  leads  to  dynamic  pulmonary  hyperinflation  during  HFOV.  The  expiratory 
flow-volume  relationship  which  leads  to  this  asymmetry  during  oscillatory  flow 
may  differ  substantially  from  that  which  prevails  during  steady  forced  maximal 
expiration.  This  difference  stems  from  an  as  yet  uncertain  mechanism,  but  may 
provide  insight  into  the  time  course  of  dynamic  airway  narrowing  or  into  the 
distribution  of  oscillatory  flow  when  it  is  further  studied. 
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TABLE  1 

Component  Values  Incorporated  into  HFOV  Models 

Values  common  to  all  3  models 

Rinsp  20  cm  E^O* sec/ liter 

Calv  0.051iter/cm  H2O 

Caw  0. 0008liter/cm  H2O 

Lbf  25  cm  sec2/ liter 

Values  pertinent  to  individual  models 

Fixed  asymmetric  resistance  model 
Rexp  40  cm  H2O/ liter 

Variable  asymmetric  resistance  model 

Rexp(Palv)  400/Palv  cm  H2O  sec/ liter  (when  Palv  is  expressed  in  cm  H2O) 

Flow  limitation  model 

Vcrit  0  cm  H20 

Rflowlim  1. 25*10^  cm  I^O’ sec/ liter 

PNP  transistor  forward  beta  1000 

These  parameters  result  in  a  maximum  expiratory  flow  rate  equal 
to  Palv*. 08  liter/ sec  (when  Palv  is  expressed  in  cm  H2O) 
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FIGURE  LEGENDS 

Figure  1.  A  schematic  representation  of  the  lung  model  used  in 
these  calculations. 

Figure  2.  Left:  Distribution  of  cross-sectional  area  as  a 
function  of  distance  from  the  upstream  (£-0)  to  downstream  (£-1) 
end.  Right:  Distribution  of  number  of  airways  located  a  distance 
£  from  the  upstream  end  of  the  airway  system. 

Figure  3.  The  dependence  of  the  loss  coefficient  CL  on  Reynolds 
number .  a 

Figure  4.  (a)  The  locus  of  points,  at  each  of  several  different 
lung  volumes,  at  which  a  transition  from  S<1  to  S>1  can  be 
achieved.  These  curves  are  obtained  by  setting  the  numerator  of 
equations  (16)  and  (17)  to  zero  with  S-l.  (b)  A  schematic 
representation  showing  how  an  elastic  jump  is  'fit'  into  the 
solution  for  Speed  index  as  a  function  of  distance.  See  text  for 
complete  description. 

Figure  5.  (a)  Dimensional  flow  rate  as  a  function  of  lung  volume 
for  several  different  combinations  of  the  parameters  of  the 
model.  Also  shown  are  several  sets  of  experimental  data  obtained 
on  humans  and  dogs.  (b)  The  location  of  the  transition  point 
(FLS)  corresponding  to  the  cases  shown  in  (a).  In  all  cases,  the 
transition  point  moves  upstream  as  lung  volume  decreases. 

Figure  6.  Plots  of  (a)  area  ratio,  (b)  speed  index,  (c)  bronchial 
pressure  minus  alveolar  pressure,  and  (d)  bronchial  pressure 
minus  atmospheric  pressure,  showing  the  distributions  with 
distance  for  several  different  levels  of  respiratory  effort. 

Figure  7.  A  collection  of  iso-volume-pressure-flow  curves  for 
different  lung  volumes  generated  by  the  model.  Also  shown  are 
the  respiratory  limits  (due  to  muscular  limitations)  as  measured 
by  Hyatt  et  al  (1966). 

Figure  8.  (a)  Dimensional  flow  rate  as  a  function  of  lung  volume 
for  several  different  gas  mixtures,  (b)  The  location  of  the  FLS 
corresponding  to  the  cases  shown  in  (a). 

Figure  9.  The  level  of  driving  pressure  (plueral  pressure) 
required  to  attain  choked  flow,  plotted  as  a  function  of  lung 
volume,  for  each  of  the  three  gas  combinations  shown  in  Figure  8. 


Figure  10.  The  distribution  of  external  pressure  introduced  to 
the  model  to  simulate  the  effect  of  a  reduction  * n  external 
pressure  for  extra-thoracic  airways. 


Figure  11.  The  effect  of  the  external  pressure  distribution  shown 
in  Figure  10  on  the  locus  of  transition  points  predicted  by  the 
model . 


Figure  12.  The  effect  of  non-uniform  external  pressure  (Figure 
10)  on  maximal  flow  rate  (a)  and  the  location  of  the  FLS  (b) . 

Figure  13-  Three  electrical  circuits  used  to  model  the  mechanical  behavior  of 
the  respiratory  system  during  HFOV.  Circuits  are  all  similar  except  for  the 
expiratory  pathways,  enclosed  by  the  interrupted  line,  (a)  Fixed  asymmetric 
resistance;  (b)  Variable  asymmetric  resistance;  (c)  Expiratory  flow  limitation. 
See  text  for  model  descriptions;  component  values  are  provided  in  Table  1. 

Figure  14.  Dependence  of  dynamic  hyperinflation  during  HFOV  on  frequency  and 
mean  airway  opening  pressure  (mean  Pao)  in  three  patients.  Ordinate  is  change 
of  lung  volume  above  FEC  with  mean  airway  opening  pressure  equal  to  atmosphere. 
Circles  with  broken  lines  denote  higher  mean  Pao;  stars  with  solid  lines  denote 
lower  mean  Pao. 

Subject  #  123 

Higher  mean  Pao  (cm  H2O)  9  8  13 

Lower  mean  Pao  (cm  H2O)  433 

Figure  15.  Dependence  of  dynamic  hyperinflation  during  HFOV  on  tidal  volume 
and  oscillatory  frequency  at  constant  mean  airway  opening  pressure  in  two 
patients.  Ordinate  is  change  of  lung  volume  above  FRC  with  mean  airway  opening 
pressure  equal  to  atmosphere.  Circles  with  broken  lines  denote  tidal  volume 
100  ml;  stars  with  solid  lines  denote  tidal  volume  50  ml. 

Figure  16.  Dependence  of  dynamic  hyperinflation  during  HFOV  on  inspiration- 
expiration  duration  ratio  (I:E)  at  constant  tidal  volume  (100  ml)  and 
11/26/85way  opening  pressure,  from  2  Series  A  patients.  Dat 

Figure  17.  Expiratory  flow  versus  lung  volume  relationships  found  in  Series  B 
patients.  The  solid  line  represents  forced  Diaximal  expiration  from  1.5  L  above 
resting  FRC;  the  data  points  shown  are  peak  expiratory  oscillatory  flow 
obtained  during  HFOV. 

Figure  18.  Predicted  dynamic  hyperinflation  during  HFOV  from  each  of  three 
nodels  of  asyo>metric  lung  impedance,  (a)  Fixed  asymmetric  resistance;  (b) 
Variable  asymmetric  resistance;  (c)  Expiratory  flow  limitation.  Circles  with 
broken  lines  denote  mean  Pao  10  cm  H2O;  stars  with  solid  lines  denote  mean 
Pao  3  cm  H2O. 

Figure  19.  Lissajous  figures  displaying  predicted  alveolar  pressure  (Palv) 
versus  predicted  airway  opening  pressure  (Pao)  for  each  of  three  models  of 
asymmetric  lung  impedance,  (a)  Fixed  asymmetric  resistance;  (b)  Variable 
asymmetric  resistance;  (c)  Expiratory  flow  limitation.  In  each  case,  Pao  leads 
Palv  by  about  one  quarter  cycle;  frequency  was  4  Hz. 

Figure  20.  Predicted  distributions  of  tidal  volume  entering  the  airway  opening 
for  each  of  the  three  models.  (A)  Central  airway  "shunt"  volume  (VTa^v);  (a) 
Fixed  asymmetric  resistance:  (b)  VariabJe  asynmetric  resistance;  (c)  Expiratory 
flow  limitation.  Circles  with  broken  lines  denote  mean  Pao  10  cm  H2O;  stars 
with  solid  lines  denote  mean  Pao  3  cm  H2O. 
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Figure  21.  Predicted  dependence  cf  distally  delivered  oscillatory  flow  (F  x 
Viaiv)  during  HFOV.  (a)  Fixed  asymmetric  resistance;  (b)  Variable  asymmetric 
resistance;  (c)  Expiratory  flow  limitation.  Circles  with  broken  lines  denote 
mean  Pao  10  cm  H2O;  stars  with  solid  lines  denote  mean  Pao  3  cm  H2O. 
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